Competition between Triplet Superconductivity and Antiferromagnetism in Quasi 

One-Dimensional Electron Systems 

Daniel Podolsky, Ehud Altman, Timofey Rostunov, and Eugene Demler 
Department of Physics, Harvard University, Cambridge MA 02138 
(Dated: February 2, 2008) 

We investigate the competition between antiferromagnetism and triplet superconductivity in quasi 
one-dimensional electron systems. We show that the two order parameters can be unified using a 
SO(4) symmetry and demonstrate the existence of such symmetry in one dimensional Luttinger 
liquids of interacting electrons. We argue that approximate SO (4) symmetry remains valid even 
when interchain hopping is strong enough to turn the system into a strongly anisotropic Fermi 
liquid. For unitary triplet superconductors SO(4) symmetry requires a first order transition be- 
tween antiferromagnetic and superconducting phases. Analysis of thermal fluctuations shows that 
the transition between the normal and the superconducting phases is weakly first order, and the 
normal to antiferromagnet phase boundary has a tricritical point, with the transition being first 
order in the vicinity of the superconducting phase. We propose that this phase diagram explains 
coexistence regions between the superconducting and the antiferromagnetic phases, and between 
the antiferromagnetic and the normal phases observed in (TMTSF^PFe. For non-unitary triplet 
superconductors the SO(4) symmetry predicts the existence of a mixed phase of antiferromagnetism 
and superconductivity. We discuss experimental tests of the SO (4) symmetry in neutron scattering 
and tunneling experiments. 

I. INTRODUCTION 

Quasi one-dimensional compounds can display a rich variety of phases, including spin-Peierls, charge density wave, 
spin density wave, and superconducting ordersiSi&iiii. Due to the large anisotropy in their crystal structure, these 
materials are often modelled as a collection of weakly coupled Luttinger liquids. The wealth of phases seen in these 
compounds is then attributed to the intrinsic instability of one dimensional electron systems towards the formation 
of quasi long range order—. As temperature is lowered, correlations along individual chains grow, until the coupling 
between chains stabilize true long range order. In the current paper, we follow this approach to study the interplay 
between triplet superconductivity (TSC) and antiferromagnetism (AF) in quasi one-dimensional electron systems. 
The starting point of our discussion is an observation that, for weak umklapp scattering, one-dimensional Luttinger 
liquids at half- filling have SO (4) symmetry at the boundary between AF and TSC phases. Near this boundary, the 
two order parameters can be unified using SO (4) symmetry, leading to strong constrains on the topology of the phase 
diagram and on the spectrum of low energy collective excitations. 

Our analysis is motivated by quasi one-dimensional Bechgaard salts (TMTSF^X, and their sulphurated counter- 
parts (TMTTF)2X. The most well studied material from this family (TMTSF^PFg is an antiferromagnetic insulator 
at ambient pressure and becomes a superconductor at high pressurei&SilfiiAi. The symmetry of the superconducting 
order parameter in (TMTSF^PFg is not yet fully established^, but there is strong evidence that electron pairing 
is spin triplet: the superconducting T c is strongly suppressed by disorde r 13 i 14 i 15 i 16 i 17 ; critical magnetic field H c2 
in the interchain direction exceeds the paramagnetic limitiSiiS; the electron spin susceptibility, obtained from the 
Knight shift measurements, does not decrease below To2£. In another material from this family, (TMTSF^ClOit, 
superconductivity is stable at ambient pressure and also shows signatures of triplet pairing 21 ' 22-23 ' 24 . Insulator to 
superconductor transition as a function of pressure has also been found for (TMTSF) 2 AsF a 2 5 . and (TMTTF) 2 PFg2&. 

There are two aspects of the SO (4) symmetry between antiferromagnetism and triplet superconductivity that we 
address in this paper. 

Classical SO (4) Symmetry. We consider the possible emergence of the classical (static) symmetry at a finite tempera- 
ture critical point. We introduce a Ginzburg-Landau (GL) free energy to describe the interaction between the AF and 
TSC orders, and we study the effects of thermal fluctuations through a large N expansion and renormalization group 
(RG) analyses in d = 4 — e and d = 2 + e dimensions. For a unitary TSC, which we argue to describe Bechgaard salts, 
we find a first order transition between AF and TSC phases, a first order transition between AF and normal phases 
ending in a tricritical point, and a weakly first order transition between TSC and normal phases. For a non-unitary 
TSC we find a mixed phase in which AF and TSC orders are present simultaneously. We argue that the system is close 
to having an SO (4) symmetric tetracritical point, but there is a narrow line of direct first order transitions between 
the normal and the mixed phase. (For a detailed discussion of the distinction between unitary and non-unitary TSC, 
see Sec. HUH) 
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Quantum SO (4) Symmetry. We introduce a quantum SO (4) rotor model which encapsulates key features of the 
competition between AF and TSC orders. We use this model to study collective excitations in the system in various 
phases. We argue that the 0-excitation, which gives one of the generators of the SO (4) algebra, should give rise to 
a sharp resonance in spin polarized neutron scattering in the TSC phase. We further predict that in the case of a 
unitary TSC the energy of the ©-resonance should decrease to nearly zero at the phase boundary with the AF phase. 
Such mode softening is not expected generally near a first order transition and would be a unique signature of the 
enhanced symmetry at the transition point. 

Bechgaard salts belong to a class of strongly correlated electron systems displaying proximity of a superconducting 
state to some kind of magnetically ordered insulating state. Other examples include the high Tc cuprates^, heavy 
fermion superconductor a 28 ' 29 , and in most cases the superconducting (SC) order parameter is spin singlet (s or d 
wave) and the insulating state has antiferromagnetic or spin density wave order. Symmetry principles have been 
introduced to study the competition of order parameters in some of these systems. In S.C. Zhang's SO(5) theory of 
high T c superconductivity^, antiferromagnetism and d-w&ve superconductivity are treated as components of a five 
dimensional order parameter. In addition to the generators of the usual charge SO (2) and spin SO (3) symmetries, new 
7r-operators are introduced, which rotate superconductivity and antiferromagnetism into each other. A combination 
of analytical approximations and numerical results can be used to argue an approximate SO (5) theory of a class 
of two dimensional lattice models, such as the Hubbard and the t-J model2i*22. The SO(5) symmetry has also 
been used to discuss quasi two-dimensional organic k-BEDT-TTF salts^. The unification approach based on higher 
symmetries has been generalized to several other types of competing states. SO (5) and SO (8) symmetries have been 
used to classify possible many-body ground states in electronic ladders^^. SO (6) symmetry has been introduced to 
discuss competing striped phases and superconductivity in the cupratess* SO(4) symmetry has been used to combine 
s-wave superconductivity and charge density wave orders in the negative U Hubbard model^L^, as well as d-wave 
superconductivity and d-density wave phase o 39 i 40 . It was also suggested that the SO (5) algebra can be used to 
combine ferromagnetism and triplet superconductivity in quasi two-dimensional S^RuO^, although the existence 
of microscopic models with such symmetry has not been demonstrated. 

There are several reasons why Bechgaard salts, and (TMTSF^PFg in particular, are promising candidates for 
experimental observation of the emergence of high symmetry from the competition of two orders. The insulator to 
superconductor transition in these materials is tuned by pressure, so the entire phase diagram can be explored in a 
single sample. This compares favorably to the cuprate superconductors, where the AF/SC transition appears as a 
function of doping and different samples are required to investigate various regimes. Another important advantage of 
Bechgaard salts is that they may be well described by a microscopic Luttingcr liquid Hamiltonian, for which we can 
demonstrate the existence of SO (4) symmetry using standard bosonization analysis. This is in contrast to the high Tc 
cuprates, in which approximations need to be made in order to even define generators of the SO(5) symmetry22ii2i^. 
A related issue is the question of quasiparticles in the AF insulating state and in the d-wave superconducting phase. 
In the former case the quasiparticle spectrum is fully gapped while in the latter case there are nodal quasiparticles. 
It is not presently known how this difference affects a quantum SO(5) symmetry for collective bosonic degrees of 
freedom. An advantage of the SO(4) symmetry in (TMTSF^PFe is that quasiparticles are fully gapped in both the 
superconducting and the insulating phases. 

Recent neutron scattering experiments demonstrated the existence of strong AF fluctuations in a triplet supercon- 
ductor SV^-RuO^*^*^. This material is not quasi one-dimensional, but it has nested pieces of the Fermi surface (see 
e.g. Ref. l44|) . Thus, we expect that this material may also show some qualitative features of the competition between 
AF and TSC discussed in this paper. 

We note that our approach is phenomenological in nature, since we do not attempt to obtain Luttinger parameters 
starting from microscopic considerations. Instead, we observe that Bechgaard salts remain strongly anisotropic even 
close to the AF/TSC phase boundary. Hence, we argue that the Luttinger parameter should be such that individual 
Id chains should be in the vicinity of such phase transition. By starting with the Luttinger Hamiltonian, we derive the 
SO (4) symmetry as its immediate consequence. We note, however, that the Luttinger liquid physics is not a necessary 
requirement for observing SO (4) symmetry near the AF/TSC phase boundary. Several groups have argued that near 
the TSC phase of Bechgaard salts, the interchain tunneling is sufficient to suppress Luttinger liquid behavior in favor 
of a strongly anisotropic Fermi liquid^S*^. We will argue below that an approximate classical SO (4) symmetry will 
be present near the AF/TSC boundary even if the ordered phases arise from a Fermi liquid state, although we still 
rely on the assumption that interchain hopping of electrons is much smaller than intrachain hopping (this condition 
is satisfied for Bechgaard salts, see Sec. IVIIl|l . Similarly, we expect that the O resonance will be present even in 
a strongly anisotropic Fermi liquid, whose observation will verify the approximate quantum SO (4) symmetry. In 
this paper, for concreteness, we will concentrate on the case where the ordered phases emerge from Luttinger liquid 
behavior on individual chains. 

This paper is organized as follows. In section [H] we discuss the Luttinger liquid model for interacting electrons in 
one dimension. For incommensurate band filling, we show that along the transition line between the TSC and the 
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FIG. 1: Phase diagram for a one dimensional system of interacting spin-1/2 fermions [50|. Here .fCp = (2-KVf + 2g4 + g\ — 
2g2)/(2nvf + 2g4 — gi + 2g2). SDW and CDW correspond to spin and charge density wave states, SS and TS to singlet and 
triplet superconducting phases. 

SDW phases, this model has SO(3)xSO(4) symmetry. At half-filling we argue that, for weak umklapp, this symmetry 
is reduced to SO (4) symmetry. For quasi one-dimensional systems such as Bechgaard salts we argue that this SO (4) 
symmetry provides a unified description of AF and TSC orders. In section ITTT1 we discuss a general GL free energy 
for the interplay between magnetism and triplet superconductivity at finite temperatures, and present mean field 
diagrams for these orders. In section Hvl we analyze thermal fluctuations using 4 — e RG analysis and demonstrate the 
absence of stable fixed points, which could control multicritical points in the phase diagram. In sectionEJwe analyze 
the case of unitary TSC competing with AF by extending the spin SO(3) group to an SO(N) algebra and using large 
N analysis. In section IVT1 we investigate the interplay of non unitary TSC and AF using large N approach and RG 
analysis for N > 3 in 4 — e and 2 + e dimensions. We also discuss a physically relevant case of N — 3. In Section 
IVIII we introduce an effective SO (4) quantum rotor model that condenses the essential features of the competition 
between the two phases. We use this model to study collective excitations in various phases. In section IVIIII we 
discuss SO(4) symmetry in highly anisotropic Fermi liquids. In section Hxl we review experimental implications of the 
SO (4) symmetry for Bechgaard salts. Finally, in section we summarize our results. 

II. MICROSCOPIC ORIGIN OF THE SYMMETRY 
A. SO(3)xSO(4) symmetry at incommensurate filling 

Consider a one dimensional electron gas with the Hamiltonian 
H = Ti-o + Hi + H-2 + 

rks 

u — 9i t t 

rti - — / , a +.k s a -,vt a +,P+qt a -,k-qs 
n, _ 92 t t 

n.2 - — 2_^a +k+qs a^ p _ qt a^. pt a + ,ks 

H4 = j;^2 a +,k+ q s a +,p-qt a +,pta+,ks + J-^2 a -,k+ q s a -,p-qt a ~., P ta~,ks (1) 

Here ks create right/left moving electrons with momenta ikf + k and spin s, and we assume linearized dispersion of 
electrons e r ^ s — (Jt = rvfk. In the Hamiltonian Q the interaction term g\ describes backward scattering and terms 52 
and 34 describe forward scattering. For now, we assume that the system has incommensurate filling, so that umklapp 
processes are not allowed. The phase diagram for this system obtained from the renormalization group analysis has 
been discussed extensively before (see e.g. Refs. and |5(J and is shown in Fig. ^ F° r the current discussion, we 
concentrate on the region of the phase diagram near the transition line between the TSC and the SDW phases at 
K p = 1, i.e. 



.9i = 2.92 



(2) 
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We demonstrate that on this line the system has an SO(3)xSO(4) symmetry that unifies order parameters of the two 
phases. 

The total spin operators are defined as 

S «=2 E a iks a ss'> a r,ks' (3) 

r,kss' 

where af s , are the usual Pauli matrices. These operators form a spin SO (3) algebra 

[S a ,Sp]=ie al3 "<S 7 (4) 
We can also combine the charge operators for right and left movers (r = ±), 

Qr = \ \ a \ks a rM ~ ~ ) (5) 

ks 

and the operators 



®l= r J2 a UA-H ( fi ) 



to form two separate isospin SO (3) algebras 

r x = i(et + e r ) 
Jy = ^(et-e r ) 

J z — Q r 

r a ,r b '] = i5 r , r ,e abc r c (?) 

The total isospin group is therefore SO(4) isoS p in =SO(3)_rxSO(3)l. Note that, since spin and isospin operators 
commute, [S a , J£] = 0, they jointly define a closed SO(3) sp j n x SO(4) isoS p in algebra. 

The total spin, S a , and the total charge, + Q— , always commute with the Hamiltonian Q). In addition, due to 
the absence of umklapp at incommensurate filling, Q + and Q- are conserved separately. As shown in Appendix lAl 
using bosonization, when the condition J2J is satisfied, the r operators also commute with the Hamiltonian. Hence, 
the system has full SO(3)xSO(4) symmetry at the phase boundary between TSC and SDW phases. We emphasize 
that the SO(3)xSO(4) symmetry of Luttinger liquids at the SDW/TSC boundary is generic and does not require fine 
tuning of the parameters. SO(4)i SOSp i n invariance has been discussed in quasi one-dimensional systems with highly 
anisotropic spin interactions^^. The O r operators in © are reminiscent of the r\ operators introduced by C.N. 
Yang to study the Hubbard model^LS, but we will show in Sec. 1111 HI that the two sets of operators define different 
symmetry groups and apply to different systems. 

Spin density wave order away from half-filling is described by a complex vector order parameter, 



$a = J2 a +,ks <7 ss' a -> k s'> ( 8 ) 



kss' 



For quasi one-dimensional systems, the band structure restricts the orbital component of the triplet superconducting 
order to be ^(p) oc p x , where x is the direction parallel to the chains. Thus, the TSC order parameter is also described 
by a complex vector, 

$ i=iE4,h(^)«"«U ( Q ) 



'.S'.S 



The factor of —i is introduced for convenience, — 102 = ( ^ J ' ^ 6 ^° Ur vec * or 0Tl ^ eT parameters Re<5, Im$, Re*, 
and Im* can be combined into a 4x3 matrix, 

(Re*),; (Im$) x (Re*), (Im% 
P = \ (Re*) y (Im#)„ (ReS) y (Im*% | (10) 
(Re*) z (Im*) 2 (Re$) z (Im*) 2 
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Each column of P transforms independently as a vector under the action of the spin group, 

67 



[S a ,P b0 ] = ie a ^P b , 



(11) 



The action of the isospin group on P is easiest to understand in terms of the operators I a = J£ ' + J a and A a = J^—J a 
For a fixed row of P, the action of the isospin generators in the basis (ReV&Q,, lm^ a , Re$ Q , Im$ Q ) is represented by 
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Q 

—i 



—i 








i 

-i 







A, 



A,, 



A. 




(12) 



Once the fourth component is identified as the "time-like" direction, this is the (Euclidean) Lorentz group, with the 
I a acting as rotations and the A a acting as boosts. Hence, the rows of order parameter Paa transform in the vector 
representation of the SO(4); sosp i n group. 



B. SO(4) symmetry at half filling 

In Bechgaard salts (TMTSF^X, three out of every four conduction states are occupied. At quarter filling, umklapp 
processes involving interactions of four electrons are allowed. Such interactions are weak, and furthermore are irrele- 
vant in the RG sense for K p > 1/4 (we remind the readers that we are interested in the regime near the SDW/TSC 
boundary, where K p 1)54. On the other hand, due to structural dimerization in Bechgaard salts^i, a gap splits 
the conduction band into a completely filled lower band and a half-filled upper band. Hence, Bechgaard salts are 
half-filled systems. At half filling, the Hamiltonian (JTJ must be modified to include two-electron umklapp scattering 
processes, 

n, 93 t t 1 33 v-* ft no\ 

= ^l2^ a +,k+qs a +,p-qt a -,pta-,k S + -^2^ a -,k+qs a -.,p-qt a +-Pt a +Ms- (13) 

Analysis of the phase diagram of Luttinger liquids at half-filling reveals that there is still a direct transition between 
AF and TSC orders at K p = 1, although this condition now corresponds 2 to g\ — 2g2 = | .^73 1 . The umklapp term 
allows scattering of two right moving electrons into two left moving ones, and vice versa. Thus, it does not commute 
with the operator A 2 = Q + — Q_, which leads to breaking of the SO(3)xSO(4) symmetry. To understand the nature 
of this symmetry breaking, it is useful to rewrite l|13fl in the form, 

Wa = H E ( Rc *(9) ' Rel(-g) - Iml(?) • Iml(- 9 )) , (14) 

where &{q) is the SDW order parameter at center of mass momentum 2k f + q, 

$ (g) = / , a+.k s Vss'a-,k-gs'- (15) 

kss' 

Equation l|14() shows explicitly that umklapp tends to pin the phase of the SDW order parameter at either or ir, 
depending on the sign of 53. This is in agreement with the observation that period two antiferromagnetic order can 
be described by a single real Neel vector. 

We will show in Section fill Bl that, whereas the Ginzburg-Landau free energy is no longer SO(3) sp j n xSO(4)j SOS pi n 
symmetric at half-filling, to linear order in 173 it maintains an SO(4)=SO(3) sp i n xSO(3)i SOS pi n symmetry. The unbroken 
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part of the isospin group, SO(3)i SOS pm, is the diagonal subgroup of SO(3)_r,xSO(3)l, which is generated by the three 
I a operators, I x = |(©1 + O), I y = ^(©^ — @), I z = Q1 where 



Q = \ ( a +,ks a +- ks + a -M a -,ks - \ 

ks 

Qt = ( a +M a +.-ki - a -M a -,-ki) ■ ( 16 ) 



k 

Without loss of generality we consider the case g% < 0, where the order parameter for antiferromagnetism is given by 
the real part of $, 

kss' 

It is easy to verify that {N, Re*, Im*} transform as vectors under both spin and isospin SO(3) symmetries. We 
define 

/(Re*), (Im*% N x \ 
Q = (Re*) y (Tm$)y N y (18) 
\ (Rc*) 2 (Im*) z N z J 

Q transforms as a vector under both SO (3) algebras 

[I a ,Qbf>] = ie abc Q c(3 

[S a ,Q b p] = ie af)l Q bl (19) 

so it describes an order parameter that transforms as a (1,1) representation of the SO(4) algebra. Since we are mostly 
interested in applying our results to Bechgaard salts, we focus mostly on this SO (4) symmetry in Ref. |5Jj as well as 
on the remainder of this paper. 

Unlike the SO(3)xSO(4) symmetry discussed in the incommensurate case, the SO(4) symmetry at half-filling is 
not a rigorous symmetry of the system. The generators of this group do not commute with the Hamiltonian of the 
system exactly. However, the main emphasis of our work is to understand the finite temperature phase diagram of 
(TMTSF) 2 PF 6 . This is obtained from the classical GL free energy, which at the AF/TSC phase boundary has SO(4) 
symmetry if we retain umklapp processes to linear order in g$ (see discussion in Section IIII B|l . In addition, with 
regards to quantum properties, SO (4) symmetry is a good starting point to study the collective modes of the system 
when c/3 is small. The latter assumption is well justified for (TMTSF^PFg, since the observed dimerization in this 
material is less than 1%£4. For small g%, modes found assuming SO (4) symmetry will have a finite overlap with 
the actual excitations of the system. In particular, the quantum numbers of the O mode discussed in Section IVIII 
including charge two and center of mass momentum 2k /, are not affected by umklapp. These properties determine 
which experimental probes couple to O. We must keep in mind, however, that the explicit breaking of SO (4) due to 
higher order corrections in 33, and also due to interchain coupling, may give a small energy gap and finite broadening to 
O, even at the AF/TSC phase boundary. We also point out that, from the point of view of O excitations, the difference 
between SO(3)xSO(4) and SO(4) symmetries corresponds to the question whether +2kf and — 2k f excitations are 
the same (SO (4) symmetry at half filling) or different (SO(3)xSO(4) symmetry away from half filling). 

Thus far in the analysis we have ignored spin-orbit effects. Microwave absorption experiments in (TMTSF^AsFg 
measured^ the anisotropy in the exchange couplings to be 10~ 6 . This ultimately determines the preferred axes for 
the Neel vector N (along the b-axis of the crystal^) and the spin component of the TSC order * (along either the a- 
or c-axisi^). However, we do not expect such tiny anisotropy to play a significant role in determining the competition 
between AF and TSC phases. We also point out that NMR experiments in (TMTSF) 2 PF 6 find a divergence of T^ 1 
at the Neel temperature that is well-described by the 0(3) isotropic Heisenberg model 58 . Thus, even for critical 
fluctuations of the AF order parameter, spin anisotropy coming from spin-orbit coupling is unobservably small. 

Before concluding this section we point out that the isospin algebra defined by equations , (JBJ , and J7J can be also 
used to relate charge density wave order and singlet superconductivity in quasi one-dimensional electron systems^. 
This is relevant for the lower half of Fig. ^ 



III. GINZBURG-LANDAU FREE ENERGY 



The main goal of this section is to investigate consequences of the SO (4) symmetry for the true finite temperature 
phase transitions, when we need to consider three dimensional fluctuations of the order parameter. One may be 
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FIG. 2: Proposed phase diagram for weak inter-chain coupling tt,. When inter-chain coupling is present, tb 7^ 0, long range 
order at finite temperatures becomes possible. In the unitary case, the second order quantum critical point of a Luttinger liquid 
becomes a first order transition between AF and unitary TSC. As tb grows, the AF phase shrinks due to reduced nesting of 
the Fermi surface. Throughout we assume positive backscattering gi > 0. 

concerned that by introducing interchain couplings, we will immediately destroy the SO (4) symmetry. As we will 
show in section lYlIII even in the case where the interchain coupling is large enough to make the system into a highly 
anisotropic Fermi liquid, approximate SO (4) symmetry prevails in the Ginzburg-Landau (GL) free energy (see e.g. Eq. 
(|29[l . This feature of the AF/TSC GL free energy implies that our analysis of the phase diagram, based on classical 
SO (4) symmetry, is valid even when the normal state is described by a highly anisotropic Fermi liquid rather than 
a collection of weakly coupled Luttinger liquids. A review of the normal state properties of organic superconductors 
at low magnetic fields is given in Ref. l47l and low temperature transport properties have been reported recently in 
Ref. IS 

We illustrate the effects of interchain coupling in Fig. [21 As the temperature is reduced on either side of K p = 1 , 
the correlation length along a chain in the appropriate correlation function grows due to intrachain interactions. At 
some finite temperature, before this length diverges, coupling between the chains becomes relevant and a true three 
dimensional transition can take place. At the crossover between one and three dimensions, the description of the 
system in terms of a Luttinger liquid is supplanted by a GL free energy describing the interactions of the order 
parameters in three spatial dimensions. In this picture, for small enough t b the presence of a phase boundary between 
AF and TSC implies that the intrachain Hamiltonian is close to the SO(4) symmetric point K p = 1. 

An important assumption of our analysis is that pressure varies the value of K p and tunes the transition across the 
AF/TSC phase boundary. We note that measurements of K p based on optical conductivity measurements have been 
carried out at ambient pressure, i.e. deep inside the SDW phase in (TMTSF) 2 PF 6 . For instance, in Refs.lrjfjlandlodl the 
value K p = 0.23 is obtained. Ref. |m] point out that this value of K p assumes that the dominant umklapp contribution 
is due to commensurability at quarter filling. They acknowledge that if umklapp is dominated by commcnsurability at 
half-filling, these measurements then imply K p = 0.925. They also point out that a half-filled model has the drawback 
of predicting a gap energy that is too small. However, it is equally difficult to justify the assumption that quarter 
filling commensurability is dominant: for K p — 0.23, the rate of divergence in the RG of commensurability at one 
quarter is exponentially smaller than the corresponding rate for commensurability at half-filling. Thus, wc feel that 
the question of the value of K p is not yet settled. 

A. Incommensurate filling 

At incommensurate filling, a translation by one lattice constant multiplies the SDW order parameter by a complex 
phase factor, 

$ ^ e 2lk f a $ = e 27,w $, (20) 

where v is the filling fraction of the conduction band. For a completely incommensurate case, when v is an irrational 
number, the GL free energy must be SO(2) symmetric with respect to the phase of <I»S2iS. In the absence of pinning 
terms, the most general GL free energy with SO(3) S pi n xSO(2) charge xSO(2)t r anslation is 

F = i|V*| 2 + i|V$| 2 + ^-|*| 2 + ^|$| 2 + Ml (|#| 2 ) 2 + U2 (|$| 2 ) 2 
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FIG. 3: Mean-field phase diagram of and for u 2 < 0. At half filling, SDW order reduces to AF order. 
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FIG. 4: Mean-field phase diagram of and j^St f° r "2 > 0. At half filling, SDW order reduces to AF order. 

+ u 3 \$ 2 \ 2 + u 4 \$ 2 \ 2 + 2 Vl \$\ 2 \$\ 2 + 2v 2 \$ ■ $\ 2 + 2v 3 \$* ■ $\ 2 (21) 

Near the phase boundary between SDW and TSC phases, for quasi one dimensional systems the form of the free 
energy is strongly constrained by the SO(3) S pi n xSO(4) isosp ; n symmetry. We expect the properties of the system to be 
well described by the free energy, 

f = \ E vp 5q vp aa + , -J2 P- aa P- aa + y E ( p L + p L - pL - pL) 

da da a 

+ Ul E Pa a P-a a P- b pP- b p + U 2 E Pa a P-apP- ba P- b p. (22) 
aa b @ aotbfi 

This is the most general free energy with SO(3) sp j n xSO(4) isospin symmetric quartic coefficients, where we have used 
the order parameter defined in equation (|10|) to display this invariance explicitly. We follow the common assumption 
that changing the external control parameters of the system only affects the quadratic coefficients. Thus, these 
are allowed to break the symmetry and tune the phase transition. For Sr ^ 0, the symmetry is broken down to 
SO(3) S p in xSO(2) chargc xSO(2) t 

ranslation ■ 

There is an explicit duality between antiferromagnetism and triplet superconductivity under reversal of the sign 
of Sr in l|22|l . The mean field phase diagram of l|22|l depends crucially on the sign of u^. For negative U2, there is 
a tendency for all vector order parameters to point along a common axis. The order parameters in this case can 
be described by a single real vector times a complex phase, ^ — e l(p n and $ = e l9 n. This is referred to in the 
3 He literature as unitary triplet superconductivity^, and in magnetism as collinear spin density order 62 . On the 
other hand, for positive U2, all vector order parameters tend to be mutually orthogonal. The real and imaginary 
parts of the order parameters can no longer be set to be parallel, Re'fxlm^ ^ and Re<f>xlm<f> ^ 0. This is the 
non-unitary / non- collinear case. The mean-field phase diagrams for (1221 for U2 negative and positive are shown in 
Figs. El and H 

In principle, the parameters of the GL free energy (|21|l can be obtained from a microscopic Hamiltonian. In 
Appendix [B] we consider a quasi one dimensional electron systems with weak interactions and obtain a free energy as 
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in with 

= 21C(3) 

167T 2 VfT 2 

The quadratic coefhcients depend on coupling constants in the TSC and SDW channels, and on temperature. They 
are typically parameterized in terms of the pressure-dependent mean-field transition temperature 

n(T,P) = a TSC (T-T c (P)) 

r 2 (T,P) = a SDW (T-T N (P)) (24) 

The analysis of Appendix [B] shows that weakly interacting Fermi liquids favor unitary TSC and collinear SDW order, 
U2 < 0. 



B. Half- filling 

As is well known, the physics of period two antiferromagnetic order is captured by a single real Neel vector. It 
is interesting to study how this comes about from the point of view of the microscopic Luttinger Hamiltonian. As 
pointed out in Section m the Hamiltonian at half- filling includes a new contribution due to umklapp scattering, Eq. 
(|14H . This gives a correction to the GL free energy, which to linear order in g% can be written as 

AF = h ((Re$) 2 - (Im$) 2 ) (25) 

where h = j^. The new term pins the SDW, and breaks the SO(4)j sosp i n =SO(3)_RxSO(3)L symmetry down to its 
diagonal subgroup SO(3)j sosp i n . In agreement with the Feynman-Hcllman theorem, the quartic coefficients derived 
in Appendix 151 are not modified to linear order in g%. Therefore, in the linear order in g^, the free energy has 
SO(4)=SO(3) S pi n xSO(3)i S ospin symmetry. In principle, the higher order contributions of the umklapp 173 term can 
break the original SO(4)i SOSp in symmetry all the way down to SO(2) c harge, generated by the total charge Q. 

The bare value of (73 in (TMTSF^PFg is small, since it is proportional to dimerization, which in this compound is 
very wealth. Assuming that Coulomb interactions are of the order of the bandwidth, this leads to a bare value of g$ of 
about 0.01. Furthermore, it is not the bare value of 173 that enters the GL free energy, but its effective (renormalizcd) 
value at the Id to 3d crossover scale. At high temperatures one-dimensional physics is observed. As temperature is 
reduced, everywhere on the TSC side of the phase diagram, as well as on the AF/TSC phase boundary, g 3 flows to 
zero. This allows us to approach the critical region from the TSC side, along which the GL free energy is SO(3)j SO 
symmetric. Even on the AF side of the phase diagram, where (73 is relevant, the flow of 53 passes near zero before 
diverging. Therefore, near the AF/TSC phase boundary, the flow spends a lot of time near zero, and the eventual 
upturn of 53 may not be reached for realistic systems, in which the 3d coupling may cut off the Id RG How at low 
temperatures. Hence, it is reasonable to take small (73 everywhere near the AF/TSC phase boundary, and to consider 
a model with SO(3)i SO symmetry. 

In what follows, we will assume that the umklapp term favors the real part of the SDW, which becomes the Neel 
order parameter N, 

N = Re$. (26) 

From now on we assume that Im$ is sufficiently well gapped, so that it does not need to be included in the analysis 
of the competition between AF and TSC. This is justified since the pinning term (|25|l is relevant in the 3d theory, 
and thus any fixed point of the theory will be characterized by strong pinning. 

It is useful to consider the relation between our SO(4) symmetry and the SO(4) symmetry introduced by C.N. Yang 
for the Hubbard model^. The symmetry generators of Yang's SO (4) is the ij operator, 

^=E c U c - tl - (27) 

k 

k summation goes over the entire Brillouin zone. This operator should be compared to our operators defined in 
eq. © for incommensurate filling, and eq. (|16|) for half-filling. Away from half filling, the difference between the 
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two operators is obvious. has momentum equal to the nesting wave vector 2k t, in contrast to 77, which always has 
center of mass momentum tt. This allows us to have SO(4) symmetry for any electron density, in contrast to Yang's 
SO (4), which only applies at half-filling. At half-filling, however, 2k f = ir, and the only difference is the relative sign 
between the left and right moving contributions. This difference is substantial. The Neel order parameter transforms 
as a singlet under the action of n, 

r i h N]=0. (28) 

This should be contrasted with O, which rotates N into the TSC order parameter Vf, see eq. (11911 . The two SO(4) 
symmetries thus differ in the order parameters which they unify, and in the microscopic models for which they apply. 
Yang's SO (4) applies to the negative U Hubbard model, for which singlet SC and CDW are degenerate lowest energy 
states at half-filling. Our SO (4) unifies AF and TSC orders, which are not degenerate for the Hubbard model. As 
we discussed earlier, we expect these to be nearly degenerate phases for half-filled systems with small umklapp (e.g. 
quarter filled systems with small dimerization, such as (TMTSF2)PFg and with K p close to one). 

In analogy with the incommensurate case, at half-filling we expect that quasi one-dimensional systems near the 
AF/TSC phase boundary have a Ginzburg-Landau free energy with SO(4)-symmetric quartic coefficients. The sym- 
metry can be made explicit in terms of the matrix order parameter Q: 



F = 



aa aa 

-u 



1 QaaQaaQbfiQbfi + ^2 QaaQafiQbaQbfi 

actbfl aafa/3 

= \ m\ 2 +l(vNr + r -±\*\ 2 + r -iN* 

+ (fi x + y)(|#| 2 ) 2 + y |*T + 2u 1 |*| 2 7V 2 + 2u 2 \$ ■ N\ 2 + (fix + u 2 )(N 2 ) 2 (29) 

Changing temperature and some other parameter of the system (e.g. pressure in (TMTSF^PFg) allows to control T\ 
and T2- The SO (4) symmetry is recovered on the line r% — r<i. 

Equation l|29|l is a special case of the most general free energy with the SO(3) x 5*0(2) symmetry of spin and charge 
rotations^, 

F = \ I V* | 2 + ^ (VAO 2 + y |*| 2 + yA? 2 + mi(|*| 2 ) 2 + u 2 (N 2 ) 2 + u 3 |*T + 2i; 1 |*| 2 Y 2 + 2v 2 \N ■ *| 2 (30) 

Translational symmetry rules out N ■ vf* x ^ because this term has a non-zero wave vector. Similarly, \N x ^I>| 2 can 
be reduced to terms already present in (|30|l . When the conditions 

n = r 2 

U 2 -U 3 = Ui 

u 2 — 2m 3 = vi 

v 2 = 2u 3 (31) 

are satisfied, we recover SO(4) symmetry. In addition, if we supplement the conditions l|31|) by 

v 2 = (32) 

there is an even higher SO(9) symmetry, which allows rotations between any components of vectors N, Re^P, and 
Im* 

F = i|V#| 2 + i(VA) 2 + ^(|*| 2 + A 2 )+iZ(|vI'| 2 +A 2 ) 2 (33) 

At half-filling, there is no distinction between collinear and non-collinear magnetism. However, the sign of u 2 still 
determines the nature of the triplet superconductivity, as well as the topology of the mean-field phase diagrams. 
These are similar to those displayed for incommensurate filling, Figs. |3]and21 the only difference being that SDW 
order is reduced to AF order. The mean-field phase diagrams for (|29|l for u 2 negative and positive are shown in 
Fig. |3| The analysis of Appendix iBl can be easily modified to half- filling, and yields an SO (4) symmetric free energy 
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of the form l(29() with coefficients still given by (|23|l . Thus, weekly interacting Fermi liquids favor the case u-i < 0. 
Strong interactions, however, can modify the quartic coefficients (|29f) . including a possible change of sign of u 2 . In 
the subsequent discussion we consider both possibilities. It is useful to note that all experimentally known cases of 
triplet pairing between fermions, such as 3 i? o 64 i 67 , and Sr2R.iJ.O4, correspond to the unitary case. Hence, negative u 2 
appears more likely. 



IV. THERMAL FLUCTUATIONS 



We now consider the free energy ((30(1 and address how fluctuations affect the mean-field phase diagram shown 
in Figs. |3I and For instance, when the quartic coefficients do not lie exactly on the SO (4) symmetric manifold, 
we will study whether such symmetry appears as we go to longer length scales and integrate out short wave length 
fluctuations. The possibility of enhanced static symmetry at the critical point has been discussed previously for 
several solid state systems. For easy axis AF in a magnetic field, SO (3) symmetry was suggested to appear at the 
spin flop critical point 6 &Si. For systems with competing singlet superconducting and antiferromagnetic orders Zhang 
suggested a static SO(5) symmetry as the bicritical point- 0j7 °. This SO(5) symmetry has also been used to study 
the quasi- two dimensional k-BEDT-TTF salts^. Yang and Zhang introduced an SO (4) symmetry for the Hubbard 
model at half-filling which unifies singlet superconductivity with charge density wave order^i^. 

To understand the role of fluctuations in models 1)30(1 and 1)29(1 we use 4 — e renormalization group (RG) analysis. 
For subsequent discussion it is useful to extend the spin SO (3) symmetry of the equation l(3U(l to a more general 
SO(N) symmetry. This is achieved by considering vectors N and $ as iV-component vectors. The RG equations can 
be derived using the standard approach^ 

^ = 2r 1 + ^L{(N+l)u 1 + 3u 3 } + ^L{Nv 1 + 2v 2 } 
dl 1 + n 1 + r 2 



dr 2 8K d , , AK d , 

-f = 2r 2 + — jL {Nv 1 +2v 2 } + — (N + 2 ) U2 

dl 1 + 7-1 l+r 2 

du\ 

~dT 



eiti - K d {(8N + 32)iti + 327ii7i 3 + 3277,3 + ANvf + 8v!V 2 + 2v 2 } 



^ = eu 2 - K d {A(N + 8)u 2 + 8Nvj + 16v 1 v 2 + 8v 2 } 

= eu 3 ~ K d {8Nu 2 3 + 48u lU3 + 2v 2 2 } 

dvi 
~dl 



ev 1 - K d {(8N + 8)7*17;! + (AN + 8)772^1 + 16t7 3 t;i + 16t7i + 8M17J2 + 4t> 2 77 2 + Av 2 } 



= ev 2 -K d {8u 1 v 2 + 8u 2 v 2 + mu 3 v 2 +32v 1 v 2 + (AN + 8)vl} (34) 

Here dl = c?A/A, where A is a momentum cut-off, and K d — ... is a surface of a unit sphere in d = 4 — e dimension. 
For the physically relevant N = 3 equations ((34(1 have only two fixed points. One is a trivial Gaussian fixed point 

»"{1,2} = ^{1,2,3} = "{1,2} = (35) 

and the other is an SO (9) Heisenberg point 

(37V + 2)e 

T{1 ' 2} ~ 6N+16 
e 

Til = 77 2 = Vl = 



(6N + 16)K d 

7*3 = v 2 = (36) 

The Gaussian fixed point is completely unstable. The SO (9) Heisenberg point has five unstable directions (for general 
N, the Heisenberg point has SO(3N) symmetry, but it remains unstable in five directions for all N > 1). The critical 
point should have only two unstable directions: 7"{i,2} should flow away from the critical point, but all the interaction 
coefficient should flow toward the fixed point. So, neither the Gaussian nor the SO(9) Heisenberg fixed points are 
good candidates for the critical point. In Fig. [5]we show RG flows in the SO(4) symmetric plane. We find two types 
of run-away flows. When we start with u 2 positive, it continues increasing. For 7i 2 negative, the RG flow make it even 
more negative. In both cases U\ flows to negative values. 
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FIG. 5: Renormalization group flow of the SO (4) symmetric theory eq. I|29fl in d = 4 — e dimensions. The sign of v.2 does not 
change under the flow, and there are no stable fixed points. Instead, there are two types of runaway flow, corresponding to 
unitary (u2 < 0) and non- unitary («2 > 0) TSC. The two are separated by a line of SO(9) symmetric theories (v,2 = 0). 

In many cases absence of a fixed point in the RG flows implies that we do not have a multicritical point in the phase 
diagram, but instead fluctuations induce a first order phase transition. Below we discuss consequences of the run-away 
flows in equations l|34|l . We point out that two types of the run-away flows in the SO (4) symmetric manifold shown 
in Fig. [S] correspond to unitary (u 2 < 0) and non-unitary (v,2 > 0) TSC. These two cases are considered separately. 



V. FINITE TEMPERATURE ANALYSIS: UNITARY CASE 

We consider model Ij29(l with iV-component vectors and with negative u 2 . In the 4 — e expansion there are run-away 
flows even for large N. Thus, we do not find a fixed point that could give a critical point. To understand the phase 
diagram in this case we employ large ./V calculations in d = 3. In the large N expansion, all bubble diagrams are 
summed self-consistentlyZSiS. The large TV approach for unitary triplet superconductors without coupling to magnetic 
order has been discussed previously ir>Z^ in the context of 3 He. 

Let us start by analyzing the superconducting phase. In the mean-field approximation the order parameter factorizes 
as \P — e ie fi. Hence, we take the average value of the order parameter in the ordered phase to be = (0, ...,0, cr) 
and separate the longitudinal and transverse components of the fluctuating part d^f = (At + iBt,Al + iBl)- For 
the Neel order parameter we also separate N — (Nt, Nl). Effective masses for At, Bt, and Nt are given by 

f A d 3 k 1 f A d 3 k 

r A = n +4(iti + u 2 )a 2 + 4(ui + u 2 )N / j- — r-T t-= VAu x N 



tb = n + 4i2icr + 4(tti + u 2 )N / 77-^7^ \- iuiN 



o (2ir) 3 k 2 +r A J (2tt) 3 \k 2 + r B k 2 + r N 

A d 3 k 1 _ „ f A d 3 k ( 1 1 



(27r) 3 fc 2 + r B J (2tt) 3 \k 2 +r A k 2 +r N 

f A d 3 k 1 f A d 3 k ( 1 1 \ 

r N = r 2 + Au^ + 4( fil +U 2 )Nj QWf 7^- + 4u lN ^ 7^3 (j^ + pT^J (3?) 

where A is the ultraviolet (short distance) cut-off of the free energy in equation (|29|l . In writing equations (|37|) we 
used that in the large N limit, ^{1,2} ~ a ~ VN , and we neglected terms of the order of 1/N, including 

contributions from longitudinal fluctuations. A requirement of the cancellation of tadpole diagrams for Al gives the 
condition r A — 0, as one would expect from the Goldstone theorem. It is convenient to define parameter r c from the 
condition 

f A d 3 k 1 

= r c + (12tii + iu 2 )N I 7—^7^2 (38) 



{^fk 2 



If we measure r's with respect to r c 



iv> = r i - r c 
tN = r 2 - r c 



(39) 



13 




we can absorb all the cut-off dependence of equations 1|37[1 into definitions of and tjsr 

/d 3 k ( 1 1 \ 

r* = t N + Au^ + Au lN J || (L 2 -^l-y^ 1+ u 2 )Nj - (40) 

Integrals are now convergent for large fc, so upper limits of integration can be extended to infinity. Solutions to 
equations (|4*U|) correspond to extremal points of free energy as a function of u. When both t$ and tjy are large, there 
are no solutions to equations l|4UI) . This is a disordered phase. Once we decrease sufficiently, a single solution 
appears at fe^M and splits into two for t^p < t^ t M- This describes the appearance of the TSC phase as a locally stable 
state. The point t^.u 1 where two solutions merge into one and disappear, correspond to the boundary of the local 
stability region of the TSC phase. As t^p is lowered further, at a temperature one of the solutions approaches 
(7 = and then disappears. This is a spinodal point below which a disordered phase is unstable to developing TSC 
order parameter. The actual first order phase transition occurs somewhere between and i^z,. 

Fig. [fJlshows a phase diagram constructed from the arguments presented above, for both TSC and AF phases. We 
note that the mixed phase with simultaneous TSC and AF orders is only possible on the = tff line. Thus, the 
first order phase transition between two types of ordered phases remains even when we include fluctuations. The 
most interesting feature of this phase diagram is that the transition between the disordered and the antiferromagnetic 
phases becomes first order in the vicinity of the critical point. 



VI. FINITE TEMPERATURE ANALYSIS: NON-UNITARY CASE 



Mean-field calculations for the free energy (|29|) with 112 > demonstrated that the TSC phase is non-unitary and 
there is a mixed phase with both TSC and AF order (see Section IIHI and figure • Within mean field theory, the 
mixed phase terminates at a tetracritical point with 5*0(4) = SO(3) x SO(3) symmetry. The goal of this section is 
to examine how the tetracritical point is affected by thermal fluctuations. 

To this end we extend the SO (3) spin symmetry to SO (TV) and approach the problem with three different methods: 
A large N analysis, a renormalization group calculation in d = 4 — e and in one in d = 2 + e. The physical picture 
that emerges form all of these approaches is that a SO(3) x SO(N) critical point exists for sufficiently large N, but 
probably does not survive down to the physical N = 3. We argue that in this case the tetracritical point is stretched 
to a line of direct first order transition from the normal state to the mixed phase. 
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FIG. 7: Phase diagram of the model l|290 with U2 > in three dimensions in the large TV limit. The four second order lines 
meet at the tetracritical point at non-zero angles. The same phase diagram appears in the d — 4 — e analysis for TV > 33 and 
in the d — 2 + e analysis for N > 5. 

A. Large TV phase diagram in three dimensions 

We consider the model (|29|) in three dimensions for large N and with u 2 > 0. We note that the quartic terms give 
a free energy that is bounded from below for u\ + u 2 /3 > 0. Thus, in this section we will always assume that this 
condition is satisfied. In appendix Owe also discuss that when u\ + u 2 /3 becomes small, of the order of 1/TV 2 (in the 
large N limit both us are of the order of 1/iV), the phase diagram may change qualitatively. 

In the mixed phase the TSC and AF both have non-zero expectation values and are orthogonal to each other. 
Hence, in the ordered phase we can choose 

(#) = (0,..., 0,^,^,0) 

(TV) = (0,..., 0,0,0,00 (41) 

Following the discussion in Section [V] we introduce longitudinal and transverse fluctuations for all order parameters. 
It is easy to verify that the requirement of cancellation of tadpole diagrams for longitudinal components implies zero 
effective masses for the transverse components. Shifting n and r 2 as in equation 1)39(1 we obtain self-consistency 
conditions for expectation values of the order parameters 

Up + (12ui + 4u 2 )cr^, + 4uict^ = 

+ (4ui + Au 2 )a 2 N + 8ui<4 = (42) 

These equations can be easily solved and we obtain a phase diagram shown in figure [7| We observe that in this case 
the only effect of fluctuations is to shift the tetracritical point from r\ = r 2 = to r± = r 2 = r c . 



B. Renormalization Group Analysis in d = 4 — e. SO(3)xSO(N) fixed point 

As shown in Fig. [S]for N = 3, all fixed points with symmetry 50(4) ~ 50(3) x SO(3) are unstable within a 4 — e 
expansion. In contrast, as N is increased a fixed point with 5*0(3) x SO(N) symmetry appears that is fully stable with 
respect to changes in the quartic interaction parameters, including those perturbations that destroy 5*0(3) x SO(N) 
symmetry. Such fixed point exists for N > 33 and has 

e 1 - 72h{N) 
~AK d N + 7 

r 2 = -2K d {(3N + 2) Vl + (N + 6)v 2 } (43) 



vi = 
v 2 = 

n = 



where u\, u 2 and U3 are related to v\ and v 2 by the constraints (|31l) . and the function h(N) = 
(N 2 + 8N - 65 + (N + 7)\/N 2 - 34A + 49)~ 1 - 1/2 A 2 + 0(1/A 3 ) is real only for N > 32. 
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The two quadratic parameters r% and r 2 are relevant, tuning the transition on a two dimensional phase diagram. 
The RG flow equations IpHI) . linearized about the fixed point ll4^|) . yield two principal directions, (Sri,Sr 2 ) oc (1,1) 
associated with the thermal exponent At, and (SrijSrz) oc (—1,2) associated with the anisotropy exponent X g . From 
these we find the critical exponents, 

* = V=l- E (4 + ^)+0(^). (44) 

Note that the crossover exponent (p for the anisotropy is always less than one. This implies^ that the phase boundaries 
meet as straight lines at the critical point, and we find the same topology of the phase diagram as shown in Fig. 



C. Renormalization Group Analysis in d = 2 + e. SO(3)xSO(N) fixed point 

The runaway flows in the equations (|34[1 mean that the system goes to strong coupling. In this limit the magnitudes 
of the vectors Reip, limp, and N have already developed locally but the directions can still fluctuate on long length 
scales. For N = 3 one of the runaway directions of equations H34[) corresponds to u 3 , v 2 > and Ui, V\ < 0. The 
corresponding strong coupling limit can be described by a triad of vectors that are all mutually orthogonal. 

F = ~J2 { R i?i( x ) ■ ei(y) + K 2 (e 2 (x) ■ e 2 (y) + e 3 (x) ■ e 3 (y))} (45) 

Here e\, e 2 , and e 3 correspond to N, Re<f , and Im$ respectively. The free energy (|45|) has an explicit SO (2) charge 
symmetry of rotations between e 2 and e 3 . The continuum version of this model is given by 



f = /^{_L ( v^ + 5 L ( 



J.(Ve 2 y + (Ve 3 y)} (46) 

Here gi oc if r and the constraints e?i • = Sij are implied. 

Let us now discuss the phase diagram of l|45|) and (|46|l . For K 12 — > we have a fully disordered phase. For 
K% t 2 -^oowe have a fully ordered phase that is a mixture of TSC and AF. When K 2 = oo the vectors e?i and e 2 are 
ordered and there is an Ising type transition between TSC and TSC+AF phases. For Ki = oo vector N is ordered 
and there is an 0(2) transition between the AF and TSC+AF states. For K 2 = there is a Heisenberg transition 
between the disordered and the AF phases. For K 3 = there is a transition between the fully disordered and the 
TSC phases. What happens in the interior of the phase diagram, however, is not clear. 

When we apply the d — 2 = e RG analysis to the model l|46|^ ^iZL2SiZ2i we obtain the flow equations 

dgi gligl + g\92 - gj) 

dl " 9l 27r( 5l + 52 )2 

— = -eg 2 + — — - — - 47 
dl 2n(g 1 + g 2 ) 2 

The flow diagram is shown in fig. |SJ We can see that it lacks the Ising and 0(2) phase transitions. This is not 
surprising, since the d — 2 = e analysis works well only for the spin- wave excitations of order parameters with TV > 3. 

To shed some light on the phase diagram of l|45(l and (|46|l we consider the large N generalization of this model. We 
assume that all vectors e*j have N components 



F = [ d d x (-^(Ver) 2 + -^((Ve 2 ) 2 + (Ve 3 ) 2 ) + -{{e, ■ Ve 2 ) 2 + (e\ • Ve 3 ) 2 ) + -(e 2 • Ve 3 ) 2 l 
J (2ffi 2.92 g 3 54 J 



(48) 



The last two terms in equation (|48l) are generated in the RG flow, even if they are absent in the microscopic model 
(such terms are linearly independent of the first two only for N > 3). The symmetry breaking pattern of the non-linear 
model (O is 



0(A) x 0(2)/0(A - 3) x 0(2) diag . 



(49) 
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In order to express the RG equations in a simple form we introduce the variables rji, 

V2 
V3 



1 




9i 




1 




92 




1 


1 


9i 


92 


2 


2 


92 


94, 



2 

93 ' 

r?4 = — + — • (50) 



These variables arise naturally in a matrix formulation of the non- linear model Eq. I|48l) . see Appendix [D] where the 
RG calculation is outlined. To one loop order we find the RG flow 

V 2 i - vZ - 7 i 

V2V3 

V2-V3- Vi _ v±_ 

2rjim 2r) 2 

N - 3 rj - vl - ill _ V±_ 
2 771772 2ij 3 



dr]i 




1 


dl 


= m - 


27 


df]2 




1 


dl 


— erj 2 — 


27 


dm 




1 


dl 


= em - 


27 


dm 




1 


dl 


= em - 


27 



N- 2 + 
— [ N-2 + 

N-2 + 



N-3r]i , rjj 



2 4 27?? 



(51) 



The conditions g\ = §2 and (73 = <?4 define a two dimensional subspace over which the free energy 148|) has the 
enhanced symmetry SO(3) x SO(N). For arbitrary AT, the RG equations (|5l"1) have an SO(3) x SO(N) fixed point 

m = T] 2 = (N- 2- x)/e, 

m = m = x vi, (52) 



where x = (N — 2 + V TV 2 — 5./V + 5)/ (N — 1). Independent of N, this point has one stable direction and one unstable 
direction within the symmetric plane. The flow in directions perpendicular to the symmetric plane depends on the 
value of N. The case TV = 4 is special and is discussed in Appendix [D] For all other N > 5, the RG flow away 
from the SO(3) x SO(N) plane has one stable and one unstable direction. Two relevant parameters are therefore 
necessary to tune the transition, just as we found in the d = 4 — e analysis for large N, eq. I|43l) . Computing the 
critical exponents v and <j>, just as we did in the case d = 4 — e, eq. (|4"4fl . we find, 
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FIG. 9: Phase diagram for systems with competing AF and non-unitary TSC orders described by the model H45|l in three 
dimensions. The tetracritical point in the mean-field phase diagram of the GL free energy in equation 129H (see Fig. |1J is 
replaced by a line of direct first order transitions between a disordered and a mixed TSC/AF phase. 



Just as before, we find that for all finite N, 4> is less than one, leading to a phase diagram that is topologically 
equivalent to that found in large N expansion and in 4 — e RG analysis (see Fig. [7J|. What's more, to leading order in 
1/./V, expansions about the upper and lower critical dimension, eqs. (|44|) and (|53(l . lead to the same critical exponents, 

l/p = rf-2 + 0(i) 

<P = l-O(^) ( 54 ) 
This supports the fact that for large N the SO(3) x SO(N) fixed point changes adiabatically with dimension. 



D. Phase Diagram for N — 3 in three dimensions 



We employed three approaches to study classical fluctuations in systems with competing AF and non-unitary TSC: 
large TV expansion, d = 4 — e and d = 2 + e RG analyses. When N is large all three consistently predict a tetracritical 
point with enhanced SO(3)xSO(N) symmetry. In the physically relevant case N — 3 and d = 3 the situation is less 
clear. For example, expansion from the upper critical dimension (d = 4 — e) shows that such fixed point appears only 
for N > 33. Expansion from the lower critical dimension (d = 2 + e) gives an SO(3)xSO(N) fixed point in the RG flow 
for any A" > 3, but these fixed points become tetracritical points on the phase diagram only for N > 5. It is possible 
that in three dimensions even for A" = 3 there is a tetracritical SO(3)xSO(3) point. The reason why perturbative 
expansions in dimension fail to see it, is that they work well for small e and extrapolations to d = 3 should be treated 
with cautionS. Such scenario, however, would contradict the results of classical Monte Carlo simulations in Ref. 
l&oL in which the model l|45l) has been analyzed for K\ — K%. In that paper the bimodal distribution in the energy 
histogram has been interpreted as a signature of the first order transition. 

The phase diagram that we propose for the model l|45|) and for systems with competing AF and non-unitary TSC 
in general is shown in Fig. Thermal fluctuations turn a tetracritical point into a line of direct first order transitions 
between a disordered and a mixed TSC/AF phase. We expect, however, that the width of such first order line is 
small and the transition is very weakly first order. We conjecture that when approaching the transition between the 
normal and the AF/TSC mixed phase, susceptibilities for the AF and TSC order parameters start increasing as if 
dominated by the SO(3)xSO(3) tetracritical point. Only very close to the transition line the divergencies are cut-off 
due to the transition being first order. 

Finally, we note that eq. I|45|) with K\ — is among a class of closely related models that have been studied 
extensively in the context of frustrated ma g ne tismSii224Si24iSSi2SiSL22i22i22i2i. Numerical studies of these models in d = 
3 dimensions yield non-universal critical properties at the boundary of normal and TSC phases^2i2&2£i2£iSL2L22i22i2&, 
and even evidence of a first order transition^. The non-perturbative theoretical analysis of Ref. supports the latter 
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scenario, claiming that the critical point observed in non-linear sigma models in d = 2 + e disappears at d c = 2.87 
in one such model, being replaced by a weakly first order transition above d c . The exact nature of the transition 
seems to be very strongly model-dependent near d = 3, and we leave open the possibility that the transition between 
non-unitary TSC and normal phases is weakly first order. 



VII. QUANTUM SO(4) SYMMETRY 



The microscopic system that motivated our discussion is an assembly of Luttinger liquids weakly coupled in three 
dimensions. It is useful to condense this system to a simpler effective quantum model that concentrates on the low 
energy collective degrees of freedom, such as AF and TSC order parameters and rotations between them (such a 
description only applies in the vicinity of the AF/TSC phase boundary shown in Fig.l). Effective quantum models 
have been discussed previously for spin systems (see Ref[93for a review), and systems with singlet superconductivity 
competing either with charge density wave orc |eB2L22iSi2i or with antiferromagnetism22i2£i2§i2L2&. 

A simple form for such an effective model is an SO (4) quantum rotor model 

'Hi = 77— X! + 7T— J! ?i _ J Qi,a.aQj,aa 
i i (ij)aa 

+ «i 5^ Qi Qi,b{3 + ^2 Qi,aaQi,a{3Qi,baQhbp 

iaba.j3 iaba.f3 

+ SrJ2(Qlza-Ql Xa -Qi,ya) (55) 

i a 

The model is obtained by coarse-graining the original lattice so that, for a half-filled system, each site of the rotor 
model includes two (or a larger even number, as necessary to include an integer number of spin-triplet Cooper pairs) 
adjacent sites along the intrachain direction of the original lattice. By combining the electronic operators that make 
up each rotor model site, one can build three local spin and three local isospin operators Si and Ii, and an SO (4) 
tensor order parameter Qi, a a- Following a procedure similar to Ref. l97i one can show that the low energy properties 
of the system are given by the rotor commutation relations, 

[Si : a , Qj.aft] — i3ij£a{3yQj.a~y 
[Ii,a: Qj.ba] — i&ij€abcQj,cot 
[Qi,aa,Qj,bp] = (56) 

These relations are analogous to equations J3J , Q , and H19|l . In equation l|55|l the unit length constraint of the rigid 
rotor models is replaced by the interaction terms u\ and U2- For Sr negative the system favors the AF state and for 
Sr positive the TSC state is preferred. When Sr = all generators of SO(4) (I = an d S = J2i ^i) commute 

with the Hamiltonian l|55|l and the system is exactly SO (4) symmetric. 

We can use equation (|55|l to discuss excitation spectra in various phases of the system. We choose to orient the AF 
order parameter (Q za ) in the z direction so that (Q zz ) — N. Similarly we take (Q xx ) — ip to describe unitary TSC, 
and (Q xx ) — (Q yy ) = ip- for non-unitary TSC. With these choices we can linearize the equations of motion for the 
fluctuations to obtain: 

^j' bl3 = b(Qbb) Ij,a^abp(Q f}p) (57) 

= -^^2^ a p'b'(Qb'b')^2(Qj,b'0' - Q h b'p'+s) ( 58 ) 

0',b> g 

-ijr = ^^2^al3'b'(Ql3'l3')^2{Qj,b'l3' -Q jth ,p, + $)+^r{ZazP>{Qp>P')Qi,zP +Za*vQj,Vz{Qzz))- ( 59 ) 



e 

&',v g 



The above equations define a linear eigenvalue problem for the frequencies of the collective modes and for the second 
quantized operators 

&t(k) = ^c T Qc 7 (k) + ]T A%S a (h) + AlhM, (60) 

c^7 a a 
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which obey £^(k) = i[H, &£] = iuj^ (k) . The Fourier transforms of lattice operators are defined by 0(k) = 
N^ 1 / 2 J2j Oj elk Xj ■ I n the effective model, neutrons couple to the spin order parameter Q za: so that the low en- 
ergy scattering intensity of polarized neutrons is given by: 

X a {k + 2k /)W ) = \(n\Q za (k)\0)\ 2 5(uj - uj n0 ) (61) 

n 

The momentum is shifted by 2k f because a uniform Q za in (|55|l corresponds to a SDW order of momentum 2k / in 
the microscopic model. The weight associated with a particular collective mode created by 6^(k) is 

\(0\b(k)Q za (k)\0)\ 2 6(uj - = \(0\[b(k),Q za (k)}\0)\ 2 5(uj - (62) 

Here we used the fact that b(k) annihilates the ground state. The commutator can be calculated using l|56|) . once the 
operator content of the mode b(k) is determined. We note that for k close to zero, neutrons couple to S a instead of 
Q za , in which case, 

xl(k,u;) = £ \{n\S a (k)\0)\ 2 S(u - u n0 ). (63) 

However, in the following, we focus on neutron scattering near 2k /. 

The nature of collective excitations in the various phases is summarized in Figs. ^J^|and Tablc||J We now provide 
a detailed analysis of the collective mode spectrum and the associated neutron scattering intensity in each phase. A 
complementary calculation of the neutron scattering intensity of the O-excitations in the unitary TSC, based on the 
microscopic model, is given in Ref. IHEL 



A. Antiferromagnet 

In the AF phase (Q zz ) = N and all other order parameters vanish. Then Eqs. (|57I59|) decouple to four independent 
collective modes. The equations of motion for the pairs {Q zx , S y } and {Q zy , S x } yield the usual AF spin waves with 
linear dispersion reflecting broken SO (3) spin symmetry: 



WAFiS (k + 2k / ) = 7vJ^(l-7 k )«4/^- |k|. (64) 



2xi V 2 X 



Here Jk = Jz/(2xi)(l — 7k), z is the lattice coordination (z = 6 for a cubic lattice in three dimensions), jk — 
z^ 1 e ik<5 , and S are bond vectors. Although k and k + 2k/ are related by reciprocal lattice vectors in the rotor 
model (|55|l . they are not in the microscopic Hamiltonian, and the addition of 2k/ to the argument of 164|) serves as a 
mnemonic for the fact that the spin mode is centered primarily around 2k/. 

Similarly the equations of motion for {Q xz ,I y } and {Q yz ,I x } describe two massive isospin waves: 



uj AFJ (k + 2k f ) = NJ^ + ^-k 2 (65) 

These excitations correspond to rotations between the AF and the TSC states and indicate proximity of the two 
ground states. When Sr goes to zero, the mass of the isospin waves vanishes reflecting an enhanced SO (4) symmetry. 
The isospin modes in the AF phase do not couple to neutrons. 



B. Unitary triplet superconductor 

In the unitary TSC (Sr > 0, 112 < 0) we choose (Q xx ) = tp, while all other order parameters vanish. Eqs. I|57(l and 
(|58|l for the pairs {Q xy , S z } and {Q XZl S y } yield two spin wave modes reflecting the broken spin symmetry: 



w«TSC,s(k) = tp^J |k|. (66) 

The creation operators for the two spin waves involve the generators 5V, and S z respectively. Substitution into 1621) 
immediately shows that the neutron scattering weight of these modes at 2k/ vanishes. 
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Eqs. (|57|l and i|59[l for the pair {Q yx ,Iz} yields a gapless phase fluctuation mode reflecting the broken charge U(l) 
symmetry in the TSC phase: 



w„TSC, v (k) = ipJ ^ |k|. (67) 

With the inclusion of Coulomb interactions, this mode becomes massive through the Higgs phenomenon, with a 
mass of the order of the plasma frequency. The creation operator of this isospin mode involves the generator I z . 
Substitution into (|62|l shows that it does not couple to 2kf neutrons. 

Finally the equations of motion for the pair {Q zx , I y } give the massive 0-mode: 



/ J Abr 

w uT sc.e(k + 2k / ) =ip\— k 2 + (68) 

V 2 X2 X2 

The operator that creates this mode from the ground state is given by 

^(k) = 7= ; , = m = r^ (Vx2(Jk 2 /2 + 4Sr)Q zx (k) + il y (k)) (69) 
V 1 + X2{Jk z /2 + Abr) \ / 

When substituted into Eq. I|62|l it gives a neutron scattering intensity at k + 2k/ : 

| (0| [6e(k), Q za (k)} \0)\ 2 b(u; - uj uTSC ,e) = 5 ™ T+j^J¥/2 + 4M 5 ( U ~ UJnTSC ^ k + 2k ^ (70) 



As we approach the point br — with SO (4) symmetry, the gap of the G-mode vanishes as v Sr. Hence, the spectrum 
of the Hamiltonian (|55(l with u?. < is such that on both sides of the AF/TSC transition we observe mode softening. 
Mode softening at the first order transition is a property of the higher symmetry quantum critical points^^^uS. 
Exactly at the SO(4) symmetric point br = the system has gapless spin and isospin wave doublets. 

C. Non unitary triplet superconductor 

For the case v,2 > the TSC phase is non-unitary. We choose (Q xx ) = (Qyy) = ip- There is also a mixed phase 
where a non-vanishing AF order parameter (Q Z z) = N appears in addition to the non unitary TSC order, considered 
below. In the pure non-unitary TSC we find two spin wave modes {Q xz , S y } and {Q yz , S x } with linear dispersion: 



W„TSC,s(k) = |k| 

As in the unitary case, these spin waves do not couple to neutrons around 2k/. 

The non unitary TSC also supports two degenerate massive 0-modes {Q zx ,I y } and {Q zy ,I x } with dispersion 



/ J A8r 

w„Tsc,e(k + 2k/) = W— fc 2 + — (72) 

V 2X2 X2 

These correspond to rotations of the real and the imaginary parts of the TSC order parameter toward the AF. These 
excitations are created by the operators: 

bljk) = (Vx2(Jky2 + 4Sr)Q zy (k) - iJ g (k)) 

V 1 + X2{Jk z /2 + Abr) \ > 

ftt (k) = 1 === ( V / X2(JfcV2 + 45r)Q zx (k) + »J y (k)) . (73) 

V 1 + X2{Jk 1 j2 + Abr) v ' 

Substituting the 6e Q operators in l|62|l we find the neutron scattering weight near 2k/: 

|(0|[6 ei (k),Q 2Q (k)]|0)| 2 ^-^ TSC ,e) = S ay 1 - — 3 f 2/2 - — b{oj - UnTscM* + 2k/)) 
\(0\[b ey (k),Q za (k)]\0)\ 2 b(u - Lo nTSC , e ) = S ax i - X2(J f 2/2 + ASr) S(w ~ ^Tsc,e(k + 2k/)) (74) 
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FIG. 10: Collective excitations in various phases, (a) In the AF phase, the excitation spectrum consists of two massless spin 
waves and two massive isospin waves. Due to translational symmetry breaking, 2k f is a reciprocal lattice vector, and these 
modes also have non-zero weight at the dashed curves near k — 0. The spectrum in the TSC phase contains massless phase 
and spin modes, as well as massive O modes. For unitary TSC, there is only one such mode. The non-unitary TSC (as well 
as the mixed phase) contains a second, degenerate mode, represented by the dotted curve, (b) The SO(4) symmetric point 
(Sr = 0). In a unitary TSC this point corresponds to the transition between AF and TSC. It is characterized by four gapless 
(Goldstone) modes (2 isospin and 2 spin). In the non-unitary case the SO(4) symmetric point is inside the mixed AF/TSC 
phase. It supports only 3 (degenerate) Goldstone modes. This is because the order parameter has a residual SO (3) symmetry 
as described in the text. 



The phase fluctuation mode in the non unitary TSC phase differs from its counterpart in the unitary case. The Eqs. 
(|57I59|) for {Q xy , Q yxi S Z ,I Z } cannot be decoupled, giving a mode that involves both spin and isospin generators. Due 
to the residual symmetry generated by S z + I z , the mode with Q xy — Q yx , drops out of the spectrum. The remaining 
excitation follows the dispersion 



U„TScA k )=^)j J (J^ + ^j k - ( 75 ) 

As in the unitary case, Coulomb interactions make this a massive mode, with a mass of the order of the plasma 
energy. The creation operator of this mode involves the generators I z and S z . Substitution into 162[l gives vanishing 
neutron scattering intensity at 2kj. 



D. Mixed phase 

In the AF/non-unitary TSC mixed phase, the order parameters form an orthogonal triad ipi = ipx, ip% = V'Z/i N = 
Nz, i.e. (Qxx) — (Qyy) = "0 an d (Qzz) = N. It is easy to verify that the phase fluctuation mode remains unchanged, 
its dispersion given by Eq. I|75(l . However, other modes are complicated due to the fact that Eqs. Ij57l59(l couple 
the coordinates {Q XZ ,Q ZX , S y , I y } and similarly {Q yz , Qzy, S x , I x }. Solution of the eigenvalue equations yields two 
collective modes for each of the above coordinate sets. One is a massive "G" mode and the other a gapless spin 
wave- like mode: 
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FIG. 11: Spin waves and O-modes in the mixed phase, (a) Dispersions (Sr > 0) of the mode clearly show mixing between the 
0-mode and the spin wave. Each of the shown modes is doubly degenerate, (b) Neutron scattering intensity of the modes at 
wave-vector k = 2kf + O.ln as a function of the tuning parameter Sr. Note that the weight of the spin wave modes goes to 
zero at the SO(4) symmetric point Sr = 0. O-modes are strongly enhanced near Sr — 0. 



Where Xt 1 = Xi 1 + X2 > = ^ z /2(l — 7k), <j) = ip 2 — N 2 , and p = ip 2 +N 2 . To calculate the spectrum in the mixed 
phase as a function of the tuning parameter Sr, we find the values of the order parameters at a given Sr from the mean 
field theory of (|55|l . Specifically we use the result ip 2 — N 2 — Srju^- Fig. Illf a - ) gives an example of the dispersions 
obtained for a particular value of Sr within the mixed phase. The asymptotic form of the excitation energies at small 
wave vectors is given by: 

~ \l^-k 2 + ^— . (77) 
V 2x2 U2X2 

Fig. ^2 demonstrates that the exact dispersion 17(j|) deviates from these asymptotic forms already at relatively small 
wave vectors. This is due to the strong mixing between spin and isospin modes. Due to this mixing, both spin 
waves and 0-modes carry some weight in the neutron scattering intensity (see Fig. ITlT b')'). The scattering intensity 
associated with the spin wave mode vanishes in the vicinity of Sr = 0. On the other hand the intensity of the 
modes becomes dramatically enhanced. Another unique feature of the phase with mixed non-unitary TSC and AF 
order is a linear with \5r\ softening of the O-excitation gap. Compare this to the y/Sr softening in the unitary TSC 
(See also Fig. El). 

The SO (4) symmetric point Sr = needs special consideration. Here N 2 = ip 2 and the order parameter is invariant 
under the SO(3) group generated by I + S. This implies that there are only three Goldstone modes at this point. 
Indeed, a direct calculation at the SO (4) symmetric point gives three degenerate modes with dispersion: 












VX. " 


X2J 



UnSOw(k)=il>J- — + — |k| (78) 



Note that the number of Goldstone modes at the SO(4) point is different in a unitary and non-unitary TSC. In 
the unitary case the spin and isospin SO(3) symmetries are broken separately with a residual U(l) x £7(1) symmetry 
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FIG. 12: Gap of the modes softens toward 8r = reflecting the enhanced SO(4) symmetry at that point, (a) The gap 
decreases as \J\5r\ in the case U2 < (unitary TSC). (b) In the case U2 > there is a change from ^/j5rj behavior in the pure 
phases to linear decrease at smaller \Sr\, inside the mixed non-unitary TSC and AF phase. 



of the order parameter, local gauge freedom associated with each. This leads to four Goldstone modes. In the non 
unitary case, on the other hand, there is residual SO(3) symmetry of the order parameter, corresponding to I + S 
rotations as discussed above. Consequently, this system has only three Goldstone modes. 

The gapless spin waves and phase modes that we found away from the SO (4) symmetry are generic to systems that 
break spin SO(3) and charge U(l) symmetries. However, excitations, which can be thought of as massive isospin 
waves, are not. Their presence shows the proximity of AF and TSC phases and their softening at the point Sr — 
should provide a unique signature of the SO (4) symmetry of the system. 



Phase 


Order parameter 


Residual symmetry 


Goldstone (massless) modes 


Pseudo-Goldstone (massive) modes 


AF 


(<?«> 


u(i)xu(i) (sr.,/,) 


2 {Sx , Sy) 


2 (/„/,) 


uTSC 




U(l) {Sx) 


3 {Sz,Sy,I Z ) 


1 {Iy) 


nuTSC 


{Qxx) = {Qyy) 


U(l) (S,+I M ) 


3 (5*2 , Sy , S Z I z) 


2 (/„,/.) 


nTSC+AF 


{Qxx) — {Qyy) ,{Qzz) 


U(l) (S z +I z ) 


3 {S Z , Sy, S Z Iz) 


2 (/„,/.) 


unitary SO (4) 


{Qxx) 


U(l)xU(l) (Ix,Sx) 


4 (Sz,Sy,Iz,I V ) 





non unitary SO (4) 


{Qxx) — {Qyy} = {Qzz} 


SO(3) (I + S) 


3 (I-S) 






TABLE I: Symmetry breaking and collective modes. Here uTSC and nuTSC stand for unitary and non-unitary TSC respectively, 
nTSC+AF corresponds to a mixed phase of non-unitary TSC and antiferromagnetism away from the SO (4) symmetric point. 



VIII. SO (4) SYMMETRY IN A STRONGLY ANISOTROPIC FERMI LIQUID 

Thus far in the analysis we have treated the case of weakly coupled Luttinger liquids, where we showed that 
SO (4) symmetry describes the phase diagram and collective modes of a system near the AF/TSC phase transition. 
However, there is a tendency away from Luttinger behavior as pressure is increased towards the superconducting state, 
as supported by the observation of field-induced SDW phases at high magnetic fieldsSSiiSfliiSiiiSSiiS&iSi, by optical 
measurements 1 ^, and by low temperature transport experiments 48 . A review of the normal state of Bechgaard salts 
at low magnetic fields is given in Ref. ^3- I n this section, we will consider the effects of interchain hopping in the 
extreme case where it is large enough to destroy all remnants of Luttinger liquid physics, and make the system into a 
highly anisotropic Fermi liquid instead. We will see that even in this limit, despite the loss of nesting, an approximate 
SO (4) symmetry remains. 

We begin by looking at the classical SO (4) symmetry of the GL free energy. In order to investigate this symmetry, 
it is sufficient to consider the quartic GL terms. This follows from the fact that our analysis of the phase diagram 
includes explicitly anisotropy in the quadratic terms, see eq. (|29H . which is tuned to zero by pressure at the AF/TSC 
transition. As shown in Appendix^] a microscopic derivation of the GL parameters starting from a weakly interacting 
Fermi liquid leads to the following form for the quartic GL terms, 



Fi = A (2(|1/| 2 ) 2 - \$ 2 f\ + B(N 2 ) 2 + 2{C + D)\$\ 2 N 2 - 4D\$ ■ N\ 2 , 



(79) 
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where, for a perfectly nested Fermi surface, A = B = C/2 = D= 16 ^2^ T 2 satisfy the SO (4) symmetry conditions H31(l 

A = D, 
B = D, 
C/2 = D, 

and the sign of the coefficients corresponds to the unitary TSC case. In the presence of interchain coupling tb, the 
single electron spectrum becomes 

£k = — 2t a cos k a — 2% cos kb — /i. (80) 

Although Bechgaard salts are triclinic, and expression ||%UJ) applies to rectangular lattices only, it gives a good approx- 
imation to the low energy quasiparticle states of the system. Here, we take t a = 250 mcV, tb = 20 meV, as estimated 
from plasma frequency measurements 1 ^ and band structure calculations 1 ^. In addition, we take \i — \[2t a , corre- 
sponding to a quarter-filled band and a nesting vector Q = (2fc/,7r) w (7r/2,7r). We note that, in the Fermi liquid 
description, dimerization only affects very high energy quasiparticles, and we exclude it from l|80|) . 

At first glance, interchain hopping seems to have a devastating effect on the SO (4) symmetry. The nesting vector 
Q no longer connects the right and left moving Fermi surfaces exactly. Hence, while the coefficient A is insensitive 
to tb, the low temperature divergence in the coefficients B, C, and D is preempted by the loss of nesting. Instead, 
these coefficients saturate at a temperature of the order of t\jt a ~ 20 K, changing the ratio A/B from unity at 
high temperatures to about 10 at T c = 1.2 K. However, nesting strongly affects antiferromagnetism only, and not 
superconductivity. Hence, it's effects on the GL parameters grow in proportion to the number of times that each GL 
parameter multiplies N in 179fl . Thus, most of the effect can be absorbed into the normalization of the field N. While 
the fields N and ^> cannot be normalized independently in the full GL free energy, as this would change the ratio of 
gradient terms ^(|V^| 2 + (ViV) 2 ), such scaling is allowed when considering the mean-field properties of the system. 
The conditions for SO (4) symmetry at mean-field level are then 

VAB = D, 

(C + D) = 3D. (81) 

Thus, at the mean-field level, SO(4) symmetry is only broken weakly. This is illustrated in Fig. EH where the left 
and right hand sides of the conditions (|81J) are evaluated explicitly. Despite the strong variation in the values of 
the different GL coefficients, the curves shown in Fig. I13l trace similar trajectories, indicating the approximate SO(4) 
symmetry. At T c = 1.2 K, we find A = 7.0 x 10 4 , B = 7.3 x 10 3 , C = 3.1 x 10 4 , and D — 2.1 x 10 4 , leading to 
VAB/D = 1.06 and (C + D)/(3D) = 0.82 (N was rescaled by a factor of 1.76). Thus, the conditions JH3J deviate 
from exact SO (4) symmetry by less than 20% at T c — 1.2K. 

While classical SO (4) symmetry is no longer exact, the phase diagram derived in previous sections does not change 
in important ways. At the mean- field level, the 20% variation from the SO (4) conditions (|81|l can open a very narrow 
mixed phase between the TSC and AF phases. Hence, as pressure is varied from the AF phase to the TSC phase, 
there will no longer be a discontinuity in the order parameters or in the density. Instead, these quantities will show 
a smooth but very rapid variation as pressure transverses the mixed phase. Therefore, due to the narrowness of the 
mixed phase, the system will be very sensitive to disorder. For realistic systems, which have impurities and crystal 
defects, the mixed phase will segregate into inhomogeneous regions of AF and TSC, just as found in the coexistence 
phase in the strictly first order case. However, unlike the case of a first order transition, the inhomogeneous behavior 
will be apparent even when the phase diagram is tuned by experimental variables that are intrinsic. We note that the 
narrowness of the mixed phase, and the corresponding sensitivity to disorder, is a direct consequence of the proximity 
of the system to SO (4) symmetry. 

We now consider the fluctuation- induced first order transition between the AF and normal phases near the SO (4) 
symmetric point. We note that in order to alter the topology of the phase diagram, the bare GL parameters must 
differ enough from the SO (4) symmetric values to divert the RG flow near a new critical point. While symmetry can 
play an important role in an RG flow, it is difficult to conceive of a situation where reduction of symmetry would lead 
to softening of the first order transitions into second order. Thus, we expect the first order transition between AF 
and normal phases discussed earlier to still be present. Finally, we note that the case we consider in this section is 
extreme, in that we look study the system as a weakly interacting, strongly anisotropic Fermi liquid. This probably 
gives a strong overestimate of the magnitude of the breaking of classical SO(4) symmetry in real systems, which are 
likely to lie between the Fermi liquid limit and the weakly coupled Luttinger liquid limit, where SO (4) is a good 
symmetry. 

Before concluding this section, we briefly discuss the effects of interchain coupling on the quantum SO (4) symmetry. 
One can look for this symmetry by verifying the existence of the O resonance in a strongly anisotropic Fermi liquid 



25 




2 4 6 8 10 



T(K) 

FIG. 13: Quartic Ginzburg-Landau coefficients l|790 for a strongly anisotropic Fermi liquid, for the choice of parameters t a = 250 
meV, tb = 20 meV. The solid curve shows \J AB, the dashed curve D, and the dash-dotted curve (C + -D)/3. Although the 
three curves do not coincide, as would be required by SO (4) symmetry, they trace similar trajectories all the way down to the 
critical temperature T c — 1.2 K (vertical dotted line). At temperatures higher than shown, the three curves converge as nesting 
is restored. 

formalism. Inside the TSC phase, this can be done usi ng R PA type calculations, which include the AF particle-hole 
and O resonance particle-particle channels (see e.g. Ref. llOSf) . Results of these calculations will be reported elsewhere. 
The main effect of interchain hopping is to fix the transverse components of the nesting vector to Q = (2k f, tt, it). 
Thus, the resonance for quasi-one dimensional systems is a collective mode, whose quantum numbers are spin zero, 
charge two, and wave vector Q. The presence of interchain coupling also introduces broadening of the mode, and 
prevents it from softening all the way down to zero energy at the AF/TSC phase transition. Instead, we expect the 
minimum energy of the excitation to be of the order of t\jt a = 16 meV. 

IX. EXPERIMENTAL SIGNATURES OF THE SO (4) SYMMETRY 

The interplay of AF and SC in the organic material (TMTSF^PFg has been a subject of active investigationiSiS. 
There is strong experimental evidence supporting that superconducting order is spin triplet, as discussed in Section^] 
In addition, (TMTSF^PFg has a quasi one-dimensional structure, as the anisotropy of electron tunneling along the 
chains (a), in the planes (b), and perpendicular to the planes (c) is of the order of t a : % : t c = 250 : 25 : 1. Hence, 
(TMTSF)2PF6 is a good candidate for comparison with the theoretical model discussed in this paper. Following the 
discussion in Section II I II and Appendix [5] we expect the triplet order parameter in this material to be unitary. The 
phase diagram for this case was obtained in Section^ In Ref. |5^we compare this to the experimental phase diagram 
of (TMTSF)2PFgi2*ii. One consequence of having enhanced symmetry at a phase transition is the suppression of the 
critical temperature due to fluctuations of one order parameter into the other. This may contribute to the drastic 
drop in Taf as pressure is increased near the AF/TSC phase boundary in Bechgaard saltsiS. 

The first order transition between AF and TSC phases near the critical point, Fig. leads to a regime of frustrated 
phase separation, with domains of one phase inside the other. The volume fractions of each phase are governed by 
the Maxwell construction, while the size of individual domains is determined by the competition between short-range 
and long-range parts of the Coulomb interactioniS*i. If the domains are distributed randomly, the total resistance of 
the system may be found using an effective medium approximation. This implies, for example, that the system is 
superconducting when the TSC phase is beyond the percolating threshold. On the other hand, it is possible that the 
TSC domains are not distributed uniformly in the system, and are more favorable on the surface of the sample. In 
this case, the TSC regions can "short-circuit" the system even before they reach the percolation condition for the 
bulk. Transport properties consistent with this scenario of an inhomogeneous system have been reported in Ref. ITol 

An interesting direction for exploring competition between AF and TSC phases is to use magnetic field experiments 
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FIG. 14: The effect of a magnetic field on the superconducting state. For points (A) far from the AF/TSC boundary, the 
magnetic field destroys superconductivity leading to a normal state. For points (B) close to the boundary, a magnetic phase 
is stabilized instead. The double line denotes a first order transition, which expands into a AF/TSC coexistence region in the 
experimental phase diagram. Here we focus on the unitary case, but similar effects can be seen for the non-unitary case near 
the AF/TSC mixed phase. Fields of the order of 100 mT are sufficient for a significant enhancement in antiferromagnetism to 
be observed. This is in contrast with field-induced SDW phases, which require fields in excess of 5 T. 



in the superconducting state near the AF/TSC phase boundary. Magnetic field produces orbital currents that strongly 
suppress electron pairing and leads to a formation of an Abrikosov vortex lattice. Suppression of the AF order by 
Zeeman effect is much smaller. Thus, we expect magnetic fluctuations to become strongly enhanced in the mixed 
state^2*S2iii2iiiiiiiS. Since the critical field along the c-axis is H£ 2 = lOOraT— , an applied magnetic field along the c- 
axis on the order of a few mT can have a strong effect, see Fig.fn] This is in contrast with effects such as field induced 
SDW's and reentrant superconductivity, which require fields of at least 5 T for their observationi*22ii22ii2iii22ii2ii24. 
For pressures close to the AF/TSC phase boundary and for slightly larger magnetic fields there may also be a quantum 
phase transition in which long range AF order develops inside the vortex phase. We note that strong sensitivity of 
1/Ti to magnetic fields in the superconducting state of (TMTSF^PFg have been reported in Ref.l20l Here, increasing 
the magnetic field along the b-axis from 12.8 mT to 232 mT results in a large increase of 1/Tj, con sistent with the 
enhancement of antiferromagnetism that we propose. Earlier specific heat measurements in Ref. Ill3l already showed 
that when the superconducting order in (TMTSF)2C104 is suppressed by a magnetic field, the system goes into a 
semimetallic state with a suppressed quasiparticle density of states. This is consistent with developing AF order, thus 
opening a gap in the quasi particle spectrum. It may be interesting to study further the enhancement of magnetic 
order in the mixed state with neutron scatteringii2iii£iii£iiii, NMR— i^L^ii, and /iSRi2i experiments. 

In Ref.|55| (see also section lVTTl of this paper) we discuss that direct observation of the 6 mode in the superconducting 
phase should be possible through neutron scattering. The most important feature of the 9-resonance, which identifies 
it as a generator of the SO (4) symmetry, is the pressure dependence of the resonance energy inside the TSC phase. 
When the pressure is reduced and the system is brought toward the phase boundary with the AF phase, we predict 
the energy of the 0-resonance to be dramatically decreased. Mode softening is not expected generically at first order 
phase transitions and provides a unique signature of the SO (4) quantum symmetry. We note that due to interchain 
hopping, the center of mass momentum of the O excitation in quasi-one dimensional systems is (2k f, n, tt). 

Another approach to detect the O excitation involves tunneling experiments with the SSC/(TMTSF)2C104 junction 
show n in Fig. 1151 (analogous experiments in the context of 7r excitations in the high Tc cuprates are discussed in Rcf. 
Il22f) . A singlet superconductor provides a reservoir of Cooper pairs that can couple to 9 pairs in (TMTSF)2C104. 
One needs to overcome, however, the momentum mismatch between the two types of pairs. A possible approach is to 
use an intermediate layer of the quasi Id material (TMTTF^PFg. This salt is quarter filled and displays spin-Peierls 
(SP) order. The modulations of the SP order thus have a periodicity of four TMTTF sites, matching the (2kf,%, 7r) 
wave vector of (TMTSF)2C104. The small mismatch between the two wave vectors, due to differences in the lattice 
constant in these compounds, can be compensated by a parallel magnetic fieldiS^. We expect peaks in the current- 
voltage characteristics of the junction when the voltage bias compensates the energy difference between Cooper and 
O pairs 

2eV = toe (82) 

Peaks in IV should be present even above the superconducting transition temperature of (TMTSF)2C104 and only 
require the other material to be superconducting. The choice of (TMTSF)2C104 is made as this material is likely to 
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FIG. 15: Tunneling experiment for detecting the excitation in (TMTSF^CICU material. A singlet superconducting material 
with a higher transition temperature than (TMTSF)2C1C>6 provides a reservoir of Cooper pairs that can couple resonantly to 
O pairs. Momentum mismatch between the Cooper pairs in SC and O pairs in (TMTSF)2C104 is compensated by scattering 
of electrons in a layer of the SP material (TMTTF) 2 PF 6 . 

be close to the AF/TSC transition at ambient pressurei^i. This eliminates the need for pressure cells, which would 
make the experiments much more difficult. 

X. SUMMARY 

The primary purpose of this paper has been to discuss the competition of antiferromagnetism and triplet super- 
conductivity in quasi one-dimensional systems, such as Bechgaard salts (TMTSF^X. The point of departure of our 
work is the existence of enhanced symmetry, that unifies the two order parameters, in one dimensional systems of 
interacting electrons. Analysis of the Luttingcr liquid model presented in Section ITU showed that the usual charge and 
spin U(l)xSO(3) symmetry is enhanced to a higher SO(3)xSO(4) symmetry on the transition line between the two 
phases for incommensurate band filling. For half-filled systems and weak umklapp scattering, the enhanced symmetry 
group becomes SO (4). Weak coupling between chains, that enables true long range order, is expected to perturb the 
SO (4) symmetry only slightly. 

In sections IlIII IIVI IVl and IVII we studied the finite temperature phase diagram for systems with SO (4) symmetry. 
For the unitary case, a mean field analysis shows that SO (4) symmetry requires a direct first order transition between 
TSC and AF phases. In addition, fluctuations of the order parameters turn a portion of the boundary between AF and 
normal phases into a first order transition and also lead to a weakly first order transition between TSC and normal 
phases. For the non-unitary case, SO(4) symmetry requires two second order transitions between TSC and AF phases. 
We find that the system is close to having an SO (4) symmetric tetracritical point, but thermal fluctuations stretch 
this point into a short line of direct first order transitions from the normal state to the mixed state. Our results have 
direct implications for quasi one-dimensional organic superconductors from the (TMTSF^X family, which are likely 
to be unitary triplet superconductors. For example, first order transitions between the AF and the TSC phases, and 
between the Normal and the AF phases explain the AF/TSC and the AF/Normal coexistence regions found in the 
phase diagram of (TMTSF) 2 PF 6 . 

In section IVTT1 we analyze collective excitations in various phases and demonstrate that SO (4) leads to the existence 
of a new collective mode, the O excitation, which describes rotations between the AF and the TSC phases. In section 
HXI we study possible experimental tests of the SO(4) symmetry. We propose that the O excitation should be observed 
as a sharp resonance in spin polarized inelastic neutron scattering experiments in the superconducting phase. We 
predict that the energy of the peak decreases toward the first order phase transition to AF order. Such softening 
of modes is not expected in general near a first order transition and would be a unique signature of the enhanced 
symmetry at the transition point. 

We thank S. Brown, P. Chaikin, M. Dressel, B.I. Halpcrin, C. Kilic, S. Sachdev, A. Turner, D.-W. Wang, and S.C. 
Zhang for useful discussions. This work was supported by Harvard NSEC. 

APPENDIX A: DERIVATION OF SO(3)xSO(4) SYMMETRY IN LUTTINGER LIQUIDS 

In this Appendix wc demonstrate that along the line g± = 2g2, the Luttinger liquid Hamiltonian has an exact 
SO(3) S pi n xSO(4)j SOSp i n symmetry. For this, we use bosonization to write the 0± operators 10 as (r = ±) 

e t = r /^^Me^l, (Al) 
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where A(x) — i\[2 ((f> p (x) + O p (x)) and 6 p = n J x dx'H p (x'). Note that are independent of the spin fields <j> a and 
6 a . Hence, the spin sector of the bosonized Hamiltonian commutes trivially with 0t, and we need only keep track of 
the charge sector, 



" 2ttK, 



(A2) 



Whenever g\ — 2g 2 , corresponding to K p = 1, the commutator [H p , A] takes on a simple form, 

[H p ,A(x)\ = V2u p (d x <p p (x)+nU p (x)) 
= -iu p d x A(x), 



(A3) 



so that commuting H p with an arbitrary function of A(x) is equivalent to taking the derivative with respect to x. For 
example, 



= \H P , A] + -{A [H p , A] + [H p , A] A{x)) + ... 



iu p [d x A + - (Ad x A + d x AA) 



= -iu p d x e J 



(A4) 



Hence, 



H p , J dxe A ^ = -iupj 



dxd x e A{x) 



= -iu a 



e A(L) _ e A(0) 



which vanishes if periodic boundary conditions are imposed on <p(x) and 6{x). Thus, for K p = 1, 



= o, 



(A5) 



(A6) 



and the Luttinger liquid has full SO(3) sp i n xSO(4)i sosp i n symmetry, generated by 0±, the total spin operators S a , and 
the charge of left and right movers, 



(A7) 



ks 



The enlarged symmetry relies on the independent conservation of total number of right- and left-movers. This is not 
a good conservation law, for instance, in the presence of impurity scattering, dimerization, or umklapp. For general 



K p , we find 



2K n 



dx (V2u p {d x ^p(x) ~ 7rII p (aO)) fJ ^ LeTA(X) - 



(A8) 



We would like to thank Daw- Wei Wang for helpful discussions on results presented in this section. 



APPENDIX B: PARAMETERS OF THE GINZBURG-LANDAU FREE ENERGY FOR WEAK 

INTERACTIONS 



To extract parameters of the GL free energy we consider a mean-field Hamiltonian 

U = y^(e fc - y)a\ s a ks + ■ Wfca^ 8 fflL fea »(o ! cr2) 8a ' + \P* ■ ^ Wfca-fc s 'afc s (o- 2 o% 

ks k k 

+ ® ■ a l-k s s a k+k ls >o ss > + & ■ ^2 al +kfS a k _ kfS ,a ss , 



(Bl) 
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FIG. 16: Diagram for the coefficient A in the GL free energy lB2t . 



where W& = |/c|/fc gives the sign of k. Integrating out the fermions yields an effective action for the order parameter 
fields. We obtain the fourth order terms: 



Fi = A[ 2(1*1 ) — |* 



>2|2 



B(2(|$| 2 ) 2 - |$ 2 | 2 



2C*|*| 2 |$| 2 + 2D ( - |$ • *f - |$* • * 



2 l/S|2 



(B2) 



where 



dk 



G 2 (-fc,-w„)G 2 (fc, W „) 



B 



C = 



D 



-y , 

2B ^ J 2tt 

2^ G2(fc ' Un)G{ - k + 2kf ' Wn) {G(fc + 2fc/ ' Wn) + 2G{k ~ 2fc/ ' Wn)} 

--3V / —G 2 {k,u n )G{-k,-u n ){G{k + 2k f , Lj n ) + G(k-2k f ,u> n )} 

J ^ G (k,Vn)G(-k,-u n )G(k + 2kf,w n )G(-h -2kf,-uj n )w- k Wk+2k f 



(B3) 



For instance, the diagram giving the coefficient is shown in Fig. 1161 For the Luttinger liquid type model with 
linearized spectrum around k — ±/c/ we obtain 



.4 



B = C/2 = D. 



(B4) 



The relationship among coefficients <|B4|1 implies that the effective GL free energy (|B2I) is SO(3) S pj n xSO(4)j sosp i n 
symmetric, as expected from the discussion in Section^] F can thus be parameterized in the form 112211 . with u\ = 3 A 
and U2 = —2A. In the clean limit, 



.4 



7C(3) 
16ir 3 v f T 2 



where £(3) = 1.202 . . ., and Vf is the Fermi velocity. 

Note that, as was pointed out in section Till Bl to linear order in c/3, umklapp does not affect the quartic coefficients 
of the free energy. For instance, the diagrams in Fig. 1171 could contribute to the coefficient of the term |<I ) Z | 2 | , I' 2 | 2 . 
However, although they do not vanish individually, the two add up to zero. This is consistent with the Feynman- 
Hellman theorem which requires that the only corrections to the free energy to linear order in be given by the 
expectation value of the perturbation (I14|) . 



APPENDIX C: NON-UNITARY TRIPLET SUPERCONDUCTIVITY AND ANTIFERROMAGNETISM. 
EXPANSION FROM THE LOWER CRITICAL DIMENSION 

Here we outline some of the methods used in the RG calculation of the non-linear model (|48|l in d = 2 + e dimensions. 
The RG flow equations of a general non-linear model can be computed using the formalism of FriedanZ^. The fields 
in such models must satisfy constrains which force them to lie on some target space manifold M.. For instance, the 
usual non-linear a model deals with a single ./V-component vector with a constrained fixed length, and M. in this case 
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FIG. 17: Corrections to coefficient of | *3?^ | 2 1 "J/z | 2 term in GL free energy, to linear order in umklapp scattering (dashed line). 
The two diagrams add up to zero. Inspection of all such diagrams shows that the quartic coefficients of the GL free energy are 
not modified to linear order. 



is the N — 1 dimensional sphere describing the locus of possible values of such vector. A local set of coordinates 4>i 
can be introduced on a patch of ./Vf, in terms of which the free energy becomes 

F = J dxg t] {4>{x))d^\x)d^{x). (CI) 

Unlike the original fields used to define the model, the fields <f>i(x) are unconstrained; all information regarding the 
original constrains is contained in the metric gij . The metric also contains the coupling constants of the system. In 
Friedan's formalism, the RG flow is thought of as a gradual deformation of the manifold as the short length degrees 
of freedom are integrated out. The RG equations can be written in a covariant way; to one loop order, 

-jjj 1 <•.!.. />'•.■ (C2) 

where e = d — 2 and Rij is the Ricci tensor, which is determined uniquely by the metric. 

In practice, whenever the manifold is a homogeneous space G/H , as in our case, it is simplest to work directly in 
the tangent space of the manifold, seeSi for a detailed discussion. In terms of the metric on the tangent space, rj a b, 
the RG equations become 

= er]ab - Rab, (C3) 
where the Ricci tensor in the tangent space is given by 

Rab = £ fac 'fib C + £ ^-^-^ /-c "he d (C4) 
U * ^ dd 

in terms of the structure factor constants of the group G 

[T a , Tb] = fab °T C + fab ! Tl 

[Ti,T b ] = f Ib c T c . (C5) 

Generators labelled by upper case indices are elements of Lie H , while lower case indices denote generators in Lie 
G-Lie H . In applying expression (|CJ4|) . we assume that the generators have been chosen so that the structure factor 
constants (|G5(1 are antisymmetric with respect to exchange of any two indices; such a choice is always possible. 
Equation l|G4|) is written in a non-covariant way to make the dependence on the coupling constants rj a b explicit, and 
it shows the advantage of working in tangent space: the Ricci tensor is given directly in terms of the Lie algebra of G. 
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We briefly digress to discuss the Lie algebra of the group SO(N), which has N(N — l)/2 generators corresponding 
to infinitesimal rotations in the planes (m, m'}, where the indices m ^ ml run through the TV independent axes. For 
instance, 5*0(3) has three generators, T x = (y,z), T y — (z,x), and T z — (x,y). Keeping in mind that (m, ml) = 
— (ml, m) ( "a clockwise rotation in the x — y plane is a counterclockwise rotation in the y — x plane" ) , we introduce a 
graphical representation: if we draw N points on a sheet of paper, an arbitrary generator (m, ml) can be represented 
by an arrow going from point m to point m! . The structure factor constants of the Lie algebra 

[{m,m'), {n,n')\ = 8 m > n (m,n') - 5 m / n >{m,n) - 5 mn (m',n') + 5 mn >(m',ri) 

can be written as a "generalized e tensor" , 

[{m, m'), (n,n')] = e{ m , m ')( n , n >) {p,p } (p, p') , (C6) 

which has a simple interpretation in terms of the arrows described above: e vanishes unless (m, m!) , (n, n') , and (p, p') 
are the edges of a closed triangle. If they do form a closed triangle, count the number of times that the directions of 
the arrows must be flipped to turn it into an oriented triangle, i.e. one satisfying m! = n, n' = p, and p' = m. If the 
number of flips is even, then e = 1; otherwise e = — 1. With this in mind, inspection of equation (|C4|) shows that, 
for groups G based on SO(N), where f ab c oc e Q i, c , the calculation of the Ricci tensor reduces almost entirely to the 
counting of triangles. 

Armed with these tools, consider the non-linear model Q48J1. 

F = f d d x (^(Vei) 2 + ^((Ve 2 ) 2 + (Ve 3 ) 2 ) + -((<?! ■ Ve 2 ) 2 + (e\ ■ Ve 3 ) 2 ) + -(e 2 • Ve 3 ) 2 l . (C7) 
J I2gi 2g 2 g 3 g 4 J 

Model l|C7|) has the symmetry SO(N) of rotations of the TV-component vectors, and the symmetry SO(2) of internal 
rotations between e*2 and S3. Hence, the symmetry group of (|C7|) is G — SO(N) x SO(2). The order parameter is 
a triad of mutually orthogonal vectors, $ = (e*i e*2 63), and the ordered phase has residual symmetry H = SO(N — 
3) x 50(2)diag- The generators of H leave the triad $ invariant, whereas the generators in Lie G-Lie H rotate the 
triad and are in one-to-one correspondence with the spin waves of the system. 

We identify four types of spin waves, corresponding to the following classes of generators: T ai , which leave e{ 2 ,3} 
untouched but rotate e*i into one of the remaining N — 3 directions; T a2 , which leave e\ untouched, but rotate either 
£? 2 or e 3 into one of the remaining N — 3 directions; T a3 , of rotations in either the ei, e 2 plane or the e*i, e*3 plane; and 
T a4 , composed of the single generator of rotations in the e*2, S3 plane. Each class furnishes an independent irreducible 
representation under the action of the group H, leading to four different spin wave velocities, and to four different 
coupling constants, r/i . . . 774, 

0,3 $00,3 ca^ ■ 



The RG flow equations HC3f) become 
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(C8) 



The fixed SO{N) x 5*0(3) symmetric point of (|C8|) is described in the body of the text for N > 5. For N — 4, 
the fixed point is stable with respect to arbitrary perturbations away from the SO(N) x 50(3) symmetric plane 
[incidentally, in this case the fixed point has a larger symmetry than expected, 50(4) x 50(4)]. On the other hand, 
within the plane, it has one stable and one unstable direction. This suggests the RG flows and the phase diagram 
that are shown schematically in Fig. El Note that we have a whole line of direct transitions between the disordered 
and the TSC+AF phases. This whole line is controlled by the point O that has a high 50(4) x 50(4) symmetry. 
This is quite remarkable: a higher symmetry appears not at a single point but at the whole transition line. 
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FIG. 18: RG flows in the equations J^TJl for N = 4. Point O has a symmetry SO (4) x SO(4). 



APPENDIX D: NON-UNITARY TRIPLET SUPERCONDUCTIVITY AND ANTIFERROMAGNETISM. 

LARGE N ANALYSIS 

In Section IVl Al we pointed out that a free energy in (|29|l in the large N limit should be considered with care, when 
fix + M2/3 is close to zero. Here we assume that 112 > 0, so, this requires negative fix. The complications arise when 
the system goes outside the basin of attraction of the tetracritical fixed point, and the RG flows carry fix to large 
negative values. As we discuss below, this leads to a first order transition which is similar to what was suggested in 
Ref. for the normal to Ai transition in liquid 3 He. We take 

"1+^/3=^ (Dl) 

where S is positive and is of order We now extend the calculations presented in Section IVI Al to the next order in 
l/N. For all order parameters we separate expectation values and fluctuations 

\P = (St + ibr, <?ip + a L + ibi, a + la^p + ib L , ax + ib ), 

N = (NtiNcNuvn + Nl), (D2) 

We can expand equation Ij29(l to order 1/A^ 2 and obtain tadpole equations for and Nl- In addition to the 
counterterms and loops due to fluctuations of the transverse components, we need to include fluctuations of the 
longitudinal components. Note, that loops of longitudinal components may be terminated by bubble chains coming 
from fi2(oT&T)ao vertices. We also need to include diagrams that arise from U2(STbT)b\a,L vertices. Special attention 
should be paid to diagonalization of propagators, since the free energy has terms which introduce mixing between 
fluctuating components in (|D2(I . 

If we want to absorb the cut-off dependence into renormalization of quadratic coefficient (compare to equation 
||3SJ), we need to define the latter relative to 



d 3 k 1 



J 



r c > = r c - (40fix + 24u 2 ) / -j^-^ — + (4fix + 8fi 2 )^ log A, (D3) 

where 



(2ir) 3 k 2 v 1 z; 2tt 2 
fi 2 N 



(D4) 



Integrals in tadpole equations can not be calculated exactly. Hence, we expand them in two cases: j 2 3> 32u2& 2 
and j 2 <C 32u2C 2 (a 2 corresponds to a 2 N or 0% or (a 2 N + a%)/2 depending on terms in the integrals). Also, while 
solving final system of equations, expansions under conditions o\, <j 2 N , a\ <C cr 2 ^, a\ ~ o 2 N were made. To be 
concrete, we took the values fix = — 1/(4JV) + 6/ (AN) and u 2 = 3/(4JV). 

Transition from disordered phase to superconductive and antiferromagnetic phases outside the vicinity of r c i remains 
of the second order, though transition border shifts such that 

15 , 16 . , 

^* = i6^v log y (D5) 

In the vicinity of full 5*0(3) x SO(N) symmetry line — we expanded equations under conditions 

\a 2 N -al\~a 2 /N, (D6) 
resulting in first order phase transition, limited by boundaries 

t N .L + 2t*,L = (D7) 
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FIG. 19: Phase diagram for large N under conditions (IDH . Solid line is a I order phase transition, dashed - II. 



and 

C 2 

tN,M + 2i*,M = 4jv2 ^ , (D8) 

where C = (1 + 3y/3/2)/ir. 

Condition (|IZ>6|) for solution obtained, appears to be valid not only for small deviations from the line of symmetry, 
but for entire line of transition. Thus solution (|D7ID8|) is self-consistent in the entire region, resulting in phase diagram 
shown on Fig. ^3 In comparison with solution of the first order expansion, boundary of mixed phase becomes a first 
order phase transition, and there is no angle between N^S <-> N and <-» >F boundaries (which is oc 5 in first order 
expansion). Boundary of the basin of attraction of stable fixed point is determined by the validity of expansion for 
different conditions for j 2 . In our case it is N5 ~ 1. 

On the other hand, expansion under condition a 2 N ^> a% for N"i£ *-* N transition, and a 2 N -C <J% for N^> «-► ^ also 
results in self-consistent solution with phase transition boundary of a different geometry. On the boundary minor 
component drops to zero, while major almost does not change. In our opinion, this solution does not have a physical 
sense. 
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